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PROBLEM  STATUS 

This  progress  report  covers  the  work  of  the 
participants  in  the  Hypervelocity  Kill  Mechanisms 
Program.  Work  on  this  problem  is  continuing. 


Authorization 
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SUMMARY 

W.W.  Atkins  -  M.A.  Persechino 
U.S.  Naval  Research  Laboratory 


INTRODUCTION 

Progress  Report  No.  17  is  a  semiannual  technical  progress 
report  covering  the  work  of  the  participants  in  the  H>->er- 
velocity  Kill  Mechanisms  Program  for  the  period  beginning 
31  March  1964  to  30  September  1964.  Reports  covering  the  work 
completed  during  and  prior  to  this  reporting  period  are  listed 
in  Section  Y. 

The  work  of  this  program  has  involved  comprehensive  studies 
designed  to  evaluate  the  feasibility  of  defeating  the  mission 
of  an  intercontinental  ballistic  missile  by  fragment  impact  and/ 
or  by  subsequent  re-entry  heating  effects.  These  effects 
include:  direct  kill  by  impact,  extent  of  aggravation  or  in¬ 
crease  in  damage  caused  by  aerothermal  effects  on  an  R/V  during 
re-entry,  aerodynamic  instability  of  nose  cones  caused  by 
damage  to  the  heat  shield  and  structure,  impact  and  thermal 
damage  to  internal  components  and  warheads,  and  perturbations 
on  the  performance  of  ICBM  booster  vehicles.  The  HKM  Program 
is  divided  into  the  following  four  phases  of  work: 

1.  Impact  Damage.  Initially  BRL,  NRL,  AVCO  and  the 
Canadian  Armament  Research  and  Development  Establishment  were 
selected  to  study  the  effects  of  hypervelocity  impacts  on  re¬ 
entry  body  materials  and  structures.  Aerojet-General  was 
selected  to  study  the  impact  effects  on  propulsion  systems. 

The  work  of  Aerojet  has  been  completed  and  the  final  report  has 
been  distributed.  The  impact  work  performed  by  AVCO  has  also 
been  completed  and  a  final  report  was  included  in  Progress  Report 
No.  13.  The  work  completed  by  CARDE  was  reported  in  Progress 
Report  No.  11.  BRL  is  preparing  a  final  report  for  their  work 
on  impacts  into  ablative  structures.  Only  NRL  is  presently 
engaged  in  impact  work  for  the  HKM  Program. 

2.  Aero thermal.  In  the  early  stages  of  the  program, 

AVCO  performed  a  multitude  of  experiments  on  cratered  heat 
shield  materials  using  rocket  exhaust  and  plasma  jet  facilities 
in  order  to  determine  the  thermodynamic  effects  on  a  damaged 
vehicle  during  re-entry .  In  the  later  stages  of  the  program, 
punctured  vehicles  (vented  and  unvented)  were  analyzed,.  GE 
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and  AVCO  performed  analytical  and  experimental  studies  on 
coupled  and  uncoupled  flows,  jet  impingement,  jet  diffusion, 
and  the  determination  of  orifice  coefficients  for  perforated 
re-entry  vehicles.  GE  conducted  an  analytical  study  to  deter¬ 
mine  the  aerodynamic  effects  on  a  damaged  vehicle  during  re¬ 
entry  (the  aeroballistic  ranges  and  the  wind  tunnels  of  NOL 
and  AEDC  were  utilized  to  provide  experimental  data).  An 
effective  kill  mechanism  did  not  evolve  from  these  studies. 
During  the  latter  part  of  the  second  year's  effort  ARAP  was 
added  to  the  participants  in  the  aerothermal  work  and,  at 
this  time,  a  strong  fundamental  research  effort  on  internal 
heating  was  established  to  determine  a  rationale  for  coupled 
and  uncoupled  flows,  impinging  jets,  and  wall  jets.  A  flight 
test  program  employing  a  NASA  propulsion  and  recovery  system 
has  been  completed  and  the  details  of  this  program  are  described 
in  Item  21  in  tbo  list  of  reports  of  Section  Y.  These  tests 
provided  both  external  and  internal  heating  data  under  actual 
and  simulated  environmental  conditions. 

3.  Vehicle  Vulnerability.  The  vulnerability  work 
initially  conducted  to  determine  the  vulnerability  of  re-entry 
body,  warhead,  and  associated  arming  and  fuzing  components 

by  BRL  and  Pica  tinny  Arsenal  has  been  terminated.  A  final 
report  on  the  vulnerability  of  nuclear  warheads  to  aero¬ 
thermal  effects  has  been  prepared  bv  Picatinny  Arsenal  and 
distributed  (See  Item  22,  Section  Y). 

Aerojet-General,  under  the  technical  management  of 
the  Weapons  Laboratories,  Detachment  4,  RTD  ,  Eglin  AFB  has 
completed  t^e  investigations  to  determine  the  vulnerability 
to  fragment  impact  of  both  liquid  and  solid  rocket  propulsion 
systems.  An  analysis  of  the  vulnerability  of  both  the  United 
States  and  other  vehicles  is  included  in  the  Aerojet  final 
report  (See  Items  24,  25  and  26). 

4.  Intelligence .  The  intelligence  phase  of  the  work  was 
designed  to  provide  information  and  guide  lines  for  the  work 
performed  in  the  other  phases  of  the  HKM  Program.  A  report 
entitled  "Soviet  ICBM  Re-Entry  Body  Study"  has  been  prepared 
by  Raytheon.  This  report  provides  a  description  of  the 
Soviet  ICBM  based  on  early  Soviet  missile  tests  in  the  Pacific 
(See  Item  1).  Additional  intelligence  data  are  described  in 
Section  T  of  previous  HKM  Progress  Reports. 

PROGRESS 

The  work  described  below  is  a  summary  of  the  technical 
progress  in  the  remaining  phases  of  the  HKM  program  fui  the 
period  ending  30  September  1964. 
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1 .  Impact  Damage  Phase 

The  investigation  by  the  Ballistic  Research  Labora¬ 
tories  of  impact  damage  to  composite  targets  utilizing  Jet 
pellets  fired  from  an  inhibited  jet  charge  Is  completed  and 
a  final  report  describing  this  work  will  lie  included  when 
available . 

The  impact  work  conducted  at  NRL  during  this  reporting 
period  has  been  directed  primarily  toward  the  development  of 
a  theoretical  understanding  of  the  mechanisms  involved  in  hyper- 
velocity  impact  damage  to  the  ablative  structure  of  re-entry 
vehicles.  Two-dimensional  wave  analysis,  utilizing  hand  cal¬ 
culations,  was  applied  to  obtain  wave  diagrams  of  shock  trans¬ 
mission  through  the  structure  and  energy  and  momentum  transfer 
from  projectile  to  target.  Pressure  -  time  recordings  were  also 
made  in  large  scale  ablative  experiments  at  impact  angles  of 
45°  in  order  to  guide  the  theoretical  approach  into  more 
realistic  impact  situations. 

In  the  experimental  aspects  of  the  program  comparisons 
were  made  of  large  and  small  scale  impacts  of  steel  and 
aluminum  projectiles  fired  into  various  thickness  phenoli.c- 
refrasil  targets  with  steel  back-up  structures.  In  general 
it  was  found  that  gross  damage  effects  on  flat-plate  ablative 
structures  can  be  extrapolated  from  smaller  scale  experiments. 
However,  extension  of  hole-size  data  to  large  scale  impacts 
at  high  angles  of  impact  indicates  a  sizable  change  in  trends. 

Three  energy  regimes  for  producing  holes  in  ablative 
structures  have  been  analyzed.  These  regimes  are  the  low 
energy  regime  where  only  the  ablative  is  perforated;  the  tran¬ 
sition  region  where  both  the  ablative  and  the  metallic  back-up 
structure  are  perforated  but  little  or  no  energy  is  dissipated 
behind  the  target;  and  the  excess  energy  region  where  the  hole 
■jize  in  the  back-up  is  equal  to  or  greater  than  the  hole  size 
in  the  ablative.  The  power  of  the  energy-thickness  ratio  in 
the  hole  size  correlation  varies  train  about:  two- thirds  to  one- 
third  depending  upon  the  energy  regimes  considered.  The  ex¬ 
cess  energy  regime  is  of  the  greater  importance  because 
impact  damage  to  components  located  inside  a  realistic,  target: 
is  considered  lethal  only  for  the  excess  energy  perforation 
condition.  It  is  significant  to  r.atn  that  in  this  regime  tin 
power  of  the  energy-thickness  term  .is  considerably  less  than 
the  power  for  this  term  in  the  transition  region,  indicating 
that  once  the  excess  energy  region  is  reached  considerably  more 
energy  per  unit  thickness  oL  material  is  required  to  pt.  educe 
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the  same  size  hole.  It  is  only  in  this  region  that  significant 
energy  is  contained  in  the  residual  or  spall  fragments  to  cause 
serious  damage  to  re-entry  vehicle  components;  and  as  a  conse¬ 
quence,  a  smaller  percentage  of  the  initial  impact  energy  is 
transferred  to  the  ablative  structure. 

In  the  scaling  experiments  conducted  for  45°  angles 
of  impact  in  the  excess  energy  regime,  the  results  show  that 
the  power  of  the  energy- thickness  term  is  greater  than  that 
for  90°  impacts  indicating  that: 

a.  At  angles  of  obliquity  larger  holes  are  produced 
in  the  ablative  than  for  normal  impacts. 

b.  More  of  the  initial  impacL  energy  is  transferred 
to  the  ablative  structure  and  less  into  residual  or  spall 
fragments.  In  these  experiments  it  was  also  noted  that  for 
ablative  thickness  to  projectile  diameter  ratios  of  0.7  more 
energy  per  unit  thickness  was  required  to  produce  the  same 
size  hole  in  structures  with  thickness  ratios  equal  to  1.4, 
indicating  a  chin  plate  effect. 

The  results  for  90°  impacts  (excess  energy  regime) 
into  astrolite  targets  with  steel  back-up  structure  were  in 
good  agreement  with  the  hole-size  correlation  presented  by  NRL 
at  the  Sixth  Symposium  for  steel  impacts  into  various  ablative 
materials . 

2 .  Aero thermal  Phase 


A  seven  test  program  was  recently  completed  in  the 
Malta  Rocket  Exhaust  Facility  to  determine  internal  heat  trans¬ 
fer  to  re-entry  vehicle  models  with  large  perforations. 

Machined  holes  were  varied  to  simulate  hypervelocity  impact 
damage  and  these  results  were  compared  with  actual  impact  damage. 
Relatively  large  perforated  models  were  tested  in  the  high 
pressure,  high  enthalpy,  ablating  environment.  Details  of  the 
first  five  tests  were  discussed  in  Section  H  of  HKM  Semiannual 
Progress  Report,  No.  15,  for  the  period  ending  March  31,  1964. 

A  summary  report  on  all  seven  tests  has  been  published  and 
appears  as  Item  29  in  the  list  of  ARPA  149  technical  reports 
in  Section  Y. 

Details  of  the  sixth  and  seventh  tests  are  presented 
in  Section  H  of  this  report  and  will  be  briefly  summarized 
below.  The  Malta  Rocket  Exhaust  Internal  Heating  Test  Program 
was  completed  with  Tests  No.  6  and  7.  The  principal  objective 
of  these  two  tests  was  to  determine  the  heat  protection 
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characteristics  of  a  lightweight  3  lb/cu  ft  rigid  urethane 
foam  when  used  as  a  filler  material  in  a  re-entry  vehicle. 
Pre-test  and  post-test  photographs  of  the  models  are  shown 
in  Figures  4  thru  7  of  Section  H.  Tempera ture- time  histories 
are  shown  in  Figures  8  and  9,  and  plots  of  foam  recession  vs 
time  are  shown  in  Figures  10  and  12. 

The  thermal  performance  of  the  foam  was  character¬ 
ized  by  an  effective  heat  of  ablation  Q*  -  £E/aM,  where 
£E  *  total  energy  entering  perforation  cavity  and 
AM  =  weight  of  foam  removed  by  ablation.  Values  of  Q*  for 
two  different  density  foams  were  determined  from  tests  Nos. 

5,  6  and  7.  These  values  are  about  2000  Btu/ lb  for  densities 
of  3.0  lb/cu  ft  and  3000  Btu/ lb  for  7.8  lb/cu  ft. 

A  principal  achievement  of  the  internal  heating 
program  was  the  establishment  of  a  correlation  of  the  rate  of 
energy  influx  through  a  heat  shield  perforation.  The  applica¬ 
tion  of  the  correlation  to  a  high-performance  unfoamed  R/V 
indicates  that  an  internal  heating  load  of  about  200  Btu  per 
square  inch  of  perforation  area  is  accumulated  by  50,000  feet 
altitude  for  a  perforation  occurring  at  100,000  feet  altitude. 
The  final  correlation  equation  for  the  rate  of  energy  ab¬ 
sorption  by  a  re-entry  vehicle  structure  is: 

dE/d0  =  0.0105  — 

y  m  l  j 

=  ratio  of  specific  heats 
P  =  local  external  static  pressure,  psia 

U  =  local  external  velocity,  ft/sec 

A  =  perforation  area,  in^ 

J  =  778  ft- lb/ Btu 

The  equation  in  this  form  neglects  the  effect  of  per¬ 
foration  enlargement  during  re-entry,  hence  calculations  for 
energy  absorption  by  the  R/V  are  considered  conservative. 

The  above  empirical  correlation  is  in  good  agreement  with  the 
results  obtained  by  Donaldson  using  a  turbulent  shear  layer 
mixing  analysis  and  is  good  over  a  range  of  three  decades  in 
the  correlating  parameter  and  for  two  widely  different  test 
environments  as  well  as  for  irregularly  shaped  holes. 

An  initial  study  was  begun  in  August  1964  with  the 
objective  of  developing  an  understanding  of  the  heat  pro¬ 
tection  characteristics  of  large  variety  of  lightweight  Loams 
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of  densities  less  than  15  lb/cu  It.  Impact  tests  have  shown 
that  foam  fillers  placed  behind  ablative  structures  are  also 
effective  in  absorbing  the  impact  energy,  therefore  the  use 
of  internal  foam  will  serve  a  dual  purpose.  Since  no  thermal 
performance  data  for  this  class  of  material,  are  presently  in 
existence,  the  study  is  being  conducted  in  two  separate  but 
related  phases.  Phase  1  is  a  concentrated  effort  on  the  basic 
understanding  of  ablation  phenomena  associated  with  cavity 
heating  and  phase  2  is  a  study  in  the  GE  5-megawatt,  air  arc 
of  the  response  of  selected  foams  with  simulated  impact 
cavities.  The  foam  response  for  two  densities  of  urethane 
foam  will  also  be  determined  in  three  different  facilities  in 
order  to  determine  the  effects  of  widely  varying  test  en¬ 
vironments. 

A  thermos Lructura 1  kill  study  is  also  being  conducted  by 
GE  to  formulate  a  simplified  computational  model  for  determin¬ 
ing  thermos  true tura 1  kill  of  perforated  re-entry  vehicles 
(both  hardened  and  unhardened).  Seven  vehicle  configurations 
with  ballistic  coefficients  from  1700  to  4400  lb/ft^  have  been 
selected  for  investigation. 

A  major  analytical  study  of  the  characteristics  of  free 
shear  flow  was  completed  by  ARAP  during  this  reporting  period. 
Expressions  have  been  derived  for  a  free  shear  flow  which 
develops  when  a  free  stream  of  gas  moving  at  some  velocity 
passes  over  a  cavity  of  stagnant  gas.  A  shear  which  develops 
between  the  moving  and  still  gas  causes  a  mixing  region  to 
develop  downstream  of  the  initial  encounter  of  the  two  gases. 

It  is  now  possible  to  solve  problems  of  this  type  with  two 
different  gases,  one  in  the  mainstream  and  the  other  in  the 
cavity.  A  principal  part  of  this  experimental  program  has  been 
the  thorough  investigation  of  jet  impingement.  Detailed  re¬ 
sults  of  this  investigation  will  be  presented  in  two  reports. 

The  first  of  these  reports  is  complete  and  appears  as  Volume  II. 
of  this  Report.  It  covers  the  free  jet  and  normal.  J-.  jingement 
studies.  A  second  report  covering  results  for  oblieje  im¬ 
pingement  is  nearing  completion  and  will  be  published  and  dis¬ 
tributed  at  a  later  date.  Most  of  the  pertinent  results  of 
both  these  reports  have  been  reviewed  in  previous  status  reports. 

The  work  reported  in  Section  K  by  AVCO  is  in  the  area  of 
instrumentation  and  investigation  of  instabilities  in  the  high 
temperature  1.5  megawatt  air  arc.  Apparatus  constructed  for 
measurement  of  gas  enthalpy,  composition,  velocity  and  density 
of  a  free  compressible,  turbulent,  high  enthalpy  jet  are  dis¬ 
cussed.  Experimental  results  from  this  high  velocity  facility 
will  be  utilized  by  ARAP  in  evaluating  the  associated  mixing 
parameter,  K,  by  means  of  the  two-gas  mixing  analysis  des¬ 
cribed  in  Section  L. 
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SUMMARY 

The  theoretical  understanding  and  prediction  of  dam¬ 
age  both  to  the  ablative  skin  and  to  the  Interior  structure  of 
a  re-entry  vehicle  for  any  impact  situation  is  currently 
being  studied  at  NRL.  In  the  experimental  phase  of  this 
program  projectile  masses  (saboted  steel  and  aluminum) 
range  from  a  few  grams  to  about  366  grams.  Impact  ener¬ 
gies  are  as  high  as  four  megajoules.  In  general  gross  dam¬ 
age  to  flat  plate  ablative  structures  extrapolates  quite  well 
from  smaller  scale  experiments;  however,  extension  of 
impact  data  at  angles  of  obliquity  into  the  megajoule  range 
indicates  a  sizable  change  in  trends. 

A  hole-size  theory  is  given  which  compares  favorably 
with  an  empirically  determined  equation.  To  reinforce  this 
study,  wave  analysis  is  used  to  deduce  relative  distribution 
between  axial  and  radial  momentum  in  the  developing  im¬ 
pact  situation.  The  theoretical  shock  wave  components  are 
studied  in  impact  experiments  where  pressure  sensors  are 
placed  at  various  points  in  and  around  the  composite  skin. 
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INTRODUCTION 

The  U.S.  Naval  Research  Laboratory  under  sponsorship  of  the  Advanced  Research 
Projects  Agency  and  Office  of  Naval  Research  has  been  engaged  In  a  study  of  hyperve¬ 
locity  impacts  in  various  re-entry  vehicle  materials  and  structural  elements.  Informa¬ 
tion  about  the  basic  deformation  process  and  formulas  predicting  damage  under  a  wide 
range  of  Impact  conditions  have  resulted.  The  objectives  of  this  work  are  twofold.  The 
first  goal  is  to  determine  the  feasibility  of  hypervelocity  fragme  *.  Impact  as  a  kill 
mechanism,  and  the  second  goal  is  to  provide  criteria  for  hardening  or  increasing  the 
resistance  of  re-entry  vehicle  structures  to  impact  damage.  The  work  discussed  in  this 
paper  is  also  directly  applicable  to  space  vehicle  technology.  1  .uch  of  the  previous 
hypervelocity  ablative  work  at  NRL  was  conducted  using  projectile  masses  of  only  a  few 
grams.  Some  doubt  as  to  the  validity  of  extrapolation  of  these  data  to  large  mass  im¬ 
pacts  has  existed  in  view  of  different  empirical  formulas  which  were  developed  from 
small  scale  work  at  various  laboratories  (1-3).  In  order  to  obta'n  data  for  large  mass 
projectiles,  an  8-inch  hypervelocity  light-gas  gun  range  was  designed  and  installed  at  the 
U.S.  Naval  Research  Laboratory  (see  Figure  1).  A  total  of  seven  shots  in  the  megajoule 
range  have  been  made  to  date  into  spaced  structure  elements  and  into  ablative  targets. 

This  paper  will  deal  with  both  large  and  small  scale  impacts  with  particular  empha¬ 
sis  on  phenolic  refrasil  targets  and  steel  and  aluminum  projectiles. 


WAVE  ANALYSIS  FOR  HYPERVELOCITY  IMPACT 
INTO  A  COMPOSITE  STRUCTURr 

The  problem  of  normal  impact  of  a  plastic  cylinder  on  ai  ve  surface  backed  up 

by  a  steel  plate  (Figure  2)  was  chosen  as  a  first  step  in  the  tlu  >  cal  consideration  for 
ablative  structure  impacts. 

This  impact  problem  is  considered  as  a  wave  propagation  problem  in  three  dimen¬ 
sions.  If  the  projectile  has  an  axis  of  symmetry  coincident  with  the  velocity  vector  and 
normal  to  the  target  surface,  then  the  problem  can  be  ideally  reduced  to  a  two-dimensional 
case.  The  methods  outlined  in  References  4  and  5  could  then  be  applied  if  the  variation 
'  •>  stre'  line  cross  sc  '  with  time  durir-  the  in*  ict  procr  is  known  (convenient 
-.nock  wave  relations  ft.,  .nmediate  characteristic  treatment  ol  the  materials  desired 
were  not  available;  in  addition,  steel  has  unique  multiple-valued  wave  transmission  prop¬ 
erties  at  high  pressure).  As  a  first  approximation,  it  was  decided  to  use  acoustic  approx¬ 
imations  of  linear  superposition  of  amplitudes  (6)  assuming  Equations  (1)  and  (2)  to  be 
usable  approximatiuns  of  the  true  physical  situation  where  ■Pi  and  .\u,  were  determined 
from  Hugoniot  curves; 
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It  is  to  be  noted  that  Equation  (2)  is  in  vector  form  and  as  inclined  waves  develop  due  to 
finite  boundary  conditions  u  and  iij  will  consist  of  both  normal  components  t^  and 
radial  components  ur.  The  Hugonlot  relating  pressure  and  the  density  function  of  oblique 
tape  wound  phenolic  nylon  (TWPN)  and  steel  were  taken  or  computed  from  References  7 
and  8.  The  impact  model  is  shown  in  Figure  2. 


Consideration  of  Tape  Wound  Phenolic  Nylon 


u,(  KM  ,  SfC) 

Figure  3  -  Shock  pressure  vs  particle  velocity 
(Hugonlot  for  pressure  determinations) 
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Figure  4  -  Pressure  vs  fractional  change 
in  density  (Hugoniot  relationships) 


for  TWPN  the  above  constants  are: 

n'  -  3.17  (km/sec)  A  =  1020  (KB) 

b'  =  1.35  K  =  3.88. 


The  relationship  connecting  shock  pressure  with  particle  velocity  for  iron  and  TWPN 
was  plotted  as  Figure  3.  The  Hugoniot  for  iron  copied  from  Reference  8  was  plotted  in 
reversed  form  (right  to  left)  from  origin  points  at  values  of  particle  velocity  u  =  3.125 
and  6.25  km/sec.  Also  the  Hugoniot  for  TWPN  was  plotted  with  initial  pressures  of  100 
and  284  KB  at  the  origin  u~  o.  Standard  procedures  were  used  to  deduce  sound  velocity 
from  the  curves  relating  pressure  and  ..  (Figure  4)  with  the  help  of  Equation  (7).  Con¬ 
sidering  the  material  as  a  compressible  fluid  the  sound  velocity  ao  is  given  by: 


since 


(7) 


(8) 


(0) 
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From  Figure  4,  P  was  determined  and  Equation  (9)  was  solved  for  the  velocity  of 
sound  over  the  range  of  values  of  .  shown  in  Figure  5  for  both  iron  and  TWPN.  An  in¬ 
teresting'  feature  of  the  curve  for  iron  is  the  discontinuity  in  sound  velocity  at  5.9  km/ sec 
exhibited  .,1  =  0.054  which  would  be  reached  at  a  pressure  of  about  130  KB  in  iron. 

When  •  is  0.21,  or  the  iron  pressure  is  about  380  KB,  a  sound  velocity  of  5.9  km/sec  is 
again  attained.  The  shaded  region  in  which  sound  velocity  falls  below  that  in  unpressur¬ 
ized  iron  can  lead  to  interesting  consequences  in  wave  propagation  as  reported  by  other 
observers  (11).  In  the  study  of  wave  transmission  through  the  steel  back-up  it  was  felt 
that  little  physical  departure  from  reality  would  be  experienced  if  the  boundary  rarefac¬ 
tion  was  assumed  to  be  a  single  shock  traveling  at  the  local  speed  of  sound  through  the 
plate.  The  plate  is  relatively  thin  compared  to  the  model  projectile  and  ablative  covering 
and  furthermore  the  spalling  criterion  which  would  result  from  a  more  exact  study  is  not 
the  objective  of  this  analysis. 

Stepwise  wave  solution  of  the  impact  model  shown  in  Figure  2  is  now  undertaken. 

For  impacts  of  like  materials  u  rl  2Tu0  3.125  km/sec.  Referring  to  Figure  3,  it  is 
observed  that  a  shock  pressure  of  284  KB  is  necessary  for  a  particle  velocity  of  3.125 
km/  sec  to  be  developed  in  stationary  TWPN.  From  Equation  (10)  (conservation  of  mass) 


Figure  5  -  Dilatational  wave  velocity  vs  density  increase  for  Fe  and  TWPN 
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<un-U)  (u-U)  ,  1.737.  (10) 

Hence,  in  the  fixed  coordinate  system,  the  shock  velocity  in  the  projectile  is  u  =  1.125 
km/sec,  and  uTWN  =  6.25  km/sec  +  1.125  km/sec  =  7.375. 

Figure  6  shows  the  shock  situation  which  has  developed  0.8  ,<sec  after  impact.  The 
wavy  line  A  in  this  figure  is  the  contact  surface  separating  the  projeclile  from  the  im¬ 
pacted  target  ablative  material.  A  compression  shock  is  traveling  back  in  the  projectile 
as  a  plane  wave  B  (edge  effects  are  neglected)  and  into  the  ablative  material  as  a  plane 
wave  designated  by  the  straight  line  C.  Since  tensile  forces  are  negligible,  a  rarefaction 
wave  (assumed  to  be  a  rarefaction  shock  to  simplify  calculation)  originates  at  the  surface 
and  proceeds  inward  radially  as  shown  by  E.  The  radial  rarefaction  originating  at  the 
junction  of  the  shock  B  and  the  cylindrical  surface  of  the  projectile  induces  a  solenoidal 
radial  expansion  motion  which  is  normal  to  the  rarefaction  wave  E  and  may  be  at  any 
angle  depending  upon  the  angle  of  the  wave  E.  The  components  of  velocity  of  the  con¬ 
tacting  materials  perpendicular  to  the  curved  surface  of  contact  are  finite  so  there  is 
an  effective  collision  velocity  of  2.22  km/sec  which  must  be  restrained  (after  initial  re¬ 
duction  of  pressure  to  zero)  in  order  to  avoid  overlap  of  projectile  and  ablative  mate¬ 
rials.  This  region  is  shown  as  a  dark  triangular  region  in  Figure  6.  A  transient  pres¬ 
sure  wave  E'  of  about  120  KB  is  generated  in  this  region  by  the  interference  and  main¬ 
tained  until  Bernoulli  flow,  necessary  to  maintain  boundary  pressure  over  the  flow  region, 
is  established.  Below  the  region  of  free  expansion  of  the  projectile  in  the  interference 
zone  there  exists  a  high  velocity  hydrodynamic  flow  region  where  the  projectile  and  tar¬ 
get  material  are  subjected  to  a  series  of  expansion  waves  until  Bernoulli  velocity  u0 


Figure  6  -  Early  stage  of  wave  study,  t  =  0,8  ,.sec 


SECRET 


SECRET 


NAVAL  RESEARCH  I  ABORATORY 


B7 


satisfies  the  Bernoulli  equation  <  0u,,2  2  ~  l>  is  reached  by  the  ejecting  material.  Com¬ 
putation  shows  u,,  to  be  6,79  km/sec. 

The  total  change  of  velocity  up  of  ejecta  when  the  ejecta  velocity  is  normal  to  the 
original  ablative  surface  is  6.25  -  (-6.79)  km/sec  or  13.04  km/sec.  A  sample  calcula¬ 
tion  on  mass  flow  of  ejecta  to  maintain  the  initial  280  KB  of  pressure  over  the  inter¬ 
ference  ring  follows  from  Newton's  law. 

PAr  2S0  gm  ■  0.81  cm-  10"  *  3.1416  -  0.62S  cm2 

W  ■  - —  - 4.  15  ■  10-’  gm  sec 

■ul>.  cm2  see2  13*  10s  cm  sec 


where 

Ar  =  Area  of  interference  ring  =  2  r  .r  =  0.625  cm2 
r  =  1.25  cm  (radius  of  projectile) 
c\r  =  0.25  cm. 

After  1.7  ..sec  (Figure  7)  the  approximate  average  mass  of  ejecta  =  w>  t  a  4.15  *  1.7  ■ 

10  1  =  0.71  gm  (where  the  mass  flow  rate  is  the  average  between  time  zero  and  1.7  ..sec). 
This  sample  calculation  was  carried  out  for  each  figurt  aid  gives  a  reasonable  check 
with  computed  mass  loss  from  the  computation  from  observed  momentum  augmentation. 
This  enlarging  propagation  pattern  continues  from  about  1.7  ..sec  when  the  compression 
wave  reaches  the  steel  back-up  surface  as  shown  in  Figure  7. 


Figure  V  -  Wave  study  (shock  arrival  at  back-up)  t  =  1.7  ..sec 
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A  sample  of  the  dynamic  values  1.7  ;<sec  after  impact  in  the  numbered  regions  of 
Figure  7  are  given  in  Table  I.  Initial  conditions  just  prior  to  impact  are  WD  =  14.9  gm, 
En  =  291  kj,  mo  =  93.1  gm-km/'sec.  At  the  laminate  rear  surface  one  has  an  interface 
problem  which  must  be  reduced  to  a  pure  impact  problem  to  apply  methods  of  Reference 
8.  This  condition  is  satisfied  if  one  pictures  a  statically  compressed  ablative  which  is 
then  impacted  by  the  steel  back-up  plate  moving  at  the  particle  velocity  u,  as  only  the 
frame  of  reference  is  changed  in  the  situation  considered.  This  situation  is  shown  in  the 
two  curves  for  TWPN  and  Fe  and  are  plotted  in  Figure  3  in  the  manner  for  impacts  with 
dissimilar  materials  with  TWPN  starting  at  the  precompressed  pressure  of  284  KB.  For 
the  area  marked  1'  in  Figure  7,  the  intersection  of  the  two  curves  in  Figure  3  is  at  610 
KB  which  is  the  predicted  pressure  wave  which  must  propagate  into  the  iron  at  about  9.3 
km/ sec.  It  is  noted  that  this  is  exactly  the  pressure  attained  if  steel  were  impacted 
directly  with  TWPN  at  6.25  km/sec.  (Right  re  intersection,  Figure  3).  This  is  greater 
than  the  284  KB  incident  wave,  so  an  additional  compressive  shock  G  of  326  KB  must  be 
propagated  back  into  the  ablative  surface  as  shown  in  inner  part  of  shock  front  F  in  Fig¬ 
ure  8.  For  the  area  between  E,  O  and  E'  and  above  F,  the  ablative  pressure  is  approxi¬ 
mately  zero.  The  intersection  of  the  bottom  TWPN  curve  of  Figure  3  is  at  about  200  KB 
which  is  the  amplitude  of  the  reflected  pressure  wave  around  G.  In  the  region  below  F 
bounded  by  E'  and  K  a  pressure  of  120  KB  is  added  to  raise  the  total  pressure  in  this 
area  to  320  KB.  The  reflected  component  F  traveling  in  a  direction  normal  to  the  target 
surface  at  a  shock  velocity  consistent  with  the  local  speed  of  sound  of  the  region  gener¬ 
ates  from  its  free  end  an  expanding  oblique  compression  wave  J  and  a  rarefaction  K  as¬ 
suming  the  applicability  of  Huygens  principle  (these  lines  assumed  to  be  straight  for 
simplification  of  calculations  would  not  exist  with  a  sharp  gradient  but  there  would  be 
rather  a  rarefaction  front  originating  at  the  reflection  boundary).  The  central  inner 


Figure  8  -  Wave  study  (reflected  shock  at 
projectile-ablative  interface)  t  =  i. 3  .ssec 
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counterpart  of  those  radial  moving  oblique  waves  are  the  respective  components  M  and  L 
originating  at  the  ends  of  G  which  are  the  junctions  of  the  326  and  200  KB  reflected  shock 
waves.  Oblique  wave  components  are  considered  negligible  between  a  200  and  280  KB 
junction  along  F.  The  radial  momentum  of  the  material  is  increased  by  this  oblique  wave 
action.  In  Figure  8  another  significant  wave  action  occurs.  The  amplitude  (negative)  of 
the  solcnoidal  converging  wave  E  rises  as  the  inverse  of  the  radius  (similar  to  a  conver¬ 
gent  cylindrical  wave)  (6)  so  that  a  point  near  the  center  is  reached  where  tensile  forces 
break  down  the  material  and  a  cavity  starts  to  form.  From  this  point  an  inner  boundary 
is  continually  formed,  and  the  wave  is  returned  as  a  compression  shock  whose  initial 
amplitude  must  fall  off  as  tile  inverse  of  tile  radius.  The  amplitude  of  the  advancing  wave 
D  also  falls  off  as  the  inverse  square  of  the  radius.  The  cavity  must  grow  as  the  mate¬ 
rial  gains  radial  momentum  outward. 

New  regions  which  arose  due  to  growing  intersecting  wave  patterns  were 
and  computations  were  made  of  their  volume,  mass,  kinetic  energy,  and  momentum 
(both  normal  and  radial).  Resultant  values  were  summed  for  each  wave  diagram  for 
graphical  presentation.  The  numbers  and  tabular  details  have  been  deleted  in  order  that 
wave  details  can  be  presented  with  greater  clarity. 

In  Figure  9,  3.3  sec  after  impact,  the  compression  wave  in  the  back-up  plate  has 
reached  the  back  boundary  and  was  reflected  as  a  rarefaction  wave  which  has  just  reached 
the  junction  between  the  ablative  and  the  back-up  accelerating  the  back-up  to  2.25  km/sec 
at  the  impact  center.  Shock  F  (originated  at  ablative-steel  boundary)  has  just  reached 
the  expanding  cone  of  the  projectile  surface  which  was  initially  cylindrical.  The  cavity 
in  the  TWPN  continues  to  grow  as  decomposed  gas  and  particles  fill  the  void  with  high 
pressure. 


Figure  ')  -  Wave  study  (cavity  formation)  t  •-  3.3  sec 
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In  the  final  wave  diagram,  Figure  10,  at  3.G  sec  after  impact  the  plate  moves  away 
from  tile  nearly  stationary  ablative  as  the  gas  in  the  cavity  continues  to  be  produced  and 
expands  into  the  void’s  with  high  pressure.  At  this  point,  the  wave  onsity  of  D  in  the 
ablative  has  fallen  (due  to  radial  attenuation)  to  about  103  KB  wher  iock  velocity  is 
less  than  that  in  steel  where  the  low  amplitude  shock  continues  at  about  7  km/sec.  This 
tends  to  loosen  the  steel  from  the  ablative  bonding  and  allows  further  pressure  action  by 
the  captive  gas  in  the  expanding  void.  The  dimensionless  subtotals  of  energy  and  mo¬ 
mentum  for  the  projectile,  ablative,  and  back-up  material  are  presented  in  graphic  form 
in  Figure  11  along  with  the  overall  totals.  These  values  are  plotted  in  terms  of  percent 
of  initial  energy  and  momentum  against  the  fractional  time  in  terms  of  the  time  it  takes 
the  projectile  to  advance  a  distance  of  one  diameter  using  its  initial  velocity.  The  energy 
and  momentum  of  the  projectile  material  falls  off  at  an  increasing  rate  as  the  projectile 
advances  into  the  target.  The  energy  and  momentum  in  the  ablative  increase  rapidly  at 
first  then  rise  more  gradually.  The  energy  in  the  ablative  falls  slightly  when  the  projec¬ 
tile  has  advanced  0.9  diameters,  =  0.9.  The  back-up  plate  momentum  rises  to  0.7  oi 
the  original  projectile  momentum  when  the  projectile  has  advanced  0.8  diameters.  About 
19f’o  of  the  original  projectile  energy  appears  in  the  back-up  material  at  ■■  0.9. 

The  total  momentum  is  about  145'o  of  Mn  which  is  about  that  to  be  expected  experi¬ 
mentally  from  impacts  into  relatively  thick,  low-strength  materials.  Total  energy  falls 
slowly  at  first  and  then  quite  steeply  to  about  65ro  of  the  original  projectile  energy  when 
the  projectile  has  advanced  0.9  diameters.  This  is  to  be  expected  since  about  10fc  of  the 
initial  kinetic  energy  has  been  carried  off  in  ejecta  which  leaves  about  25'1  of  the  energy 
in  the  ablative  and  trapped  gases  in  the  form  of  potential  energy,  which  is  still  increasing. 

In  Figure  12  consideration  is  next  given  to  the  increase  of  the  value  of  scalar  sum¬ 
mation  of  radial  momentum  in  the  projectile,  ablative,  and  back-up  material.  The  growth 


Figure  10  -  Wave  study  (back-up  ablative  sepa  rat  inn)  t  3,0  sec 
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Figure  11  -  Energy  and  momentum  vs  time  (wave  study  summations) 


in  the  value  of  this  scalar  summation  of  radial  momentum  is  most  rapid  in  the  ablative 
due  to  immediate  high-pressure  radial  propagation  of  the  initial  shock.  A  maximum  of 
62.5  gm-km/sec  of  scalar  radial  momentum  is  reached  when  the  projectile  has  advanced 
0.85  D0.  The  rise  in  scalar  radial  momentum  in  the  projectile  is  less  rapid,  reaching 
slightly  over  50  gm-km/sec  when  the  projectile  has  advanced  0.9  do  .  It  appears  that  the 
back-up  may  attain  a  maximum  of  20  gm-km/sec  sometime  later.  In  Figure  13  consid¬ 
eration  is  given  to  the  growth  of  total  scalar  radial  momentum  whose  rise  is  at  first 
gradual  with  50%  of  the  original  value  reached  when  the  projectile  has  advanced  0.5  D0. 
The  growth  in  scalar  radial  momentum  is  then  quite  rapid  with  a  maximum  of  nearly 
140%  M0  being  reached  after  advance  of  0,9  of  the  projectile.  This  is  in  agreement  with 
radial  momentum  in  spall  observations  on  moderately  thick,  low-strength  targets  (1,9). 

To  conclude  the  wave  portion  of  the  study,  it  is  noted  that,  in  spite  of  extremely  sim¬ 
plifying  assumptions  (i.e.,  rarefaction  shock,  linear  superposition,  amplitudes  propor¬ 
tioned  per  Inverse  of  radius)  the  wave  study  gives  good  results  in  predicting  momentum 
multiplication,  observed  forward  ejection  of  material  back  along  the  trajectory  at  high 
velocity,  and  growth  in  radial  momentum  leading  to  opening  of  the  target  skin.  The  growth 
of  a  large  pocket  of  heated  ablative  gas  between  the  ablative  and  back-up  which  will  later 
propel  bits  of  the  shattered  ablative  back  along  the  trajectory  is  certainly  in  accord  with 
photographic  observations  of  impac  ts  on  ablative  targets.  There  can  be  no  doubt  of  the 
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Figure  1  i  -  Increase  of  radial  momentum  with  time 


extensive  chemical  decomposition  of  the  phenolic  compounds  under  the  action  of  high 
amplitude  multiple  shocks.  The  extensive  deposition  of  condensed  carbon  following  a 
massive  ablative  shot  indicates  that  there  was  abundant  carbon  freed  from  hydrogen- 
carbon  bonding  so  there  should  also  be  released  corresponding  amounts  of  free  hydrogen. 

Seven  models  of  the  cracking  out  of  ablative  beyond  the  extreme  pressure  zone  close 
to  the  impact  point  have  been  considered.  One  model  of  a  developing  petaling  action  and 
ablative  shearout  is  shown  in  Figure  14.  The  trapped  gas  cavity  B  developed  by  wave 
processes  (also  shown  in  Figure  10)  has  continued  to  expand  and  to  eject  material  through 
the  central  impact  zone  where  small  ruptures  developed  in  both  the  ablative  and  back-up 
structure.  The  ablative  and  back-up  plate  are  pushed  apart  by  the  high  pressure  gas 
developed  in  the  immediate  impact  zone,  by  ablative  breakdown  into  gaseous  constituents 
and  by  the  energy  absorbed  from  the  high-pressure  shock.  It  is  assumed  that  though 
small  cracks  develop  in  the  accelerating  ablative,  they  interlock  so  that  gas  is  stUl  con¬ 
tained  and  forces  are  transmitted  to  maintain  an  effective  plate  area  for  the  gas  to  act 
against.  The  force  K  developed  by  gas  pressure  P  is  held  by  shear  resistance  s  of  the 
ablative  in  the  growing  gas  pressurized  ablative  disk  until  some  critical  diameter  equal 
to  D„  is  reached  at  which  cracking  and  shearout  of  the  shock  conditioned  material  occurs. 
This  situation  may  be  described  mathematically  by  Equation  (11). 
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Since  En  is  Lhe  only  source  of  energy  to  release  ablative  gas,  this  leads  to  expressions 
in  Equation  (12),  where  there  is  sufficient  excess  energy  to  cause  rupture  of  the  back-up 
in  time  to  allow  hemispherical  gas  expansion  prior  to  ablative  failure 


PV,,:  k ,  and  v(,  k  1  nj(' . 


Solving  Equations  (11)  and  (12)  and  l)n  gives  Equation  (13). 
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There  are  transition  cases  where  the  back-up  does  not  rupture  in  time  to  allow  ef¬ 
fective  hemispherical  gas  expansion  because  the  shock  is  overtaken  by  the  rarefactions 
originating  first  at  the  sides  and  then  at  the  back  of  the  projectile.  Usually  in  this  regime 
the  plate  is  given  a  strong  shock  pulse  which  is  attenuated  at  some  stage  before  the  shock 
can  make  the  entire  traverse  through  the  plate  and  then  back  as  a  rarefaction  to  enable 
the  plate  to  acquire  the  full  velocity  necessary  to  quickly  achieve  a  full  petaling  type  of 
rupture  needed  to  achieve  the  hemispherical  type  of  gas  expansion  postulated  in  the  ex¬ 
cess  energy  theory  for  ablative  hole  formation.  Deep  rounded  bowing  of  the  back-up  in¬ 
cluding  moderate  petaling  occurs  in  the  final  stages  of  this  energy-thickness  regime  for 
a  given  ablative  thickness.  It  appears  that  in  this  case  shear  and  viscous  forces  restrain 
the  high  velocity  center  and  a  high  velocity  transverse  wave  is  generated  which  moves 
out  radially  augmenting  the  impulse  from  the  ablative  wave  so  that  in  the  early  stages  a 
disk-shaped  cavity  is  opened  such  that  the  volume  can  be  approximated  by 


v  -  —  d  . 

K  a  un 


(14) 


Solving  Equations  (11),  (12)  and  (14)  as  before  gives  Equation  (15)  for  the  low  energy 
transition  regime. 


k 4  T„> 
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EXPERIMENTAL  RESULTS 

When  the  wave  model  was  conceived  the  only  complete  set  of  Hugoniot  curves  avail¬ 
able  for  ablative  materials  was  for  TWPN.  The  wave  development,  however,  is  valid 
qualitatively  for  any  ablative  up  to  Equation  (11),  except  for  metallic  projectiles  where 
there  will  be  a  set  of  reflection  interactions  which  will  raise  the  pressure  of  the  inter¬ 
posed  ablative  by  a  set  of  step  increments,  This  is  shown  experimentally  in  pressure 
recordings  in  a  later  discussion.  The  multiple  extreme-pressure  shocks  can  raise  the 
temperature  of  the  interposed  region  in  excess  of  decomposition  temperatures  for  most 
ablative  components.  Most  phenolics  under  this  condition  will  release  hydrogen  and 
gaseous  carbon  in  amounts  which  will  vary  directly  with  the  size  of  the  impact  and  also 
with  ablative  characteristics.  Choosing  astrollte  as  an  ablative  material  fixes  s  and  K, 
to  values  valid  only  for  astrollte.  In  the  preliminary  examination  concerned  with  getting 
a  gross  answer  covering  the  general  class  of  ablatives,  the  scatter  in  data  was  as  great 
as  the  variations  in  the  factors 
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so  that  for  1  ...i.  ■-  o  considerations,  the  best  lit  for  the  whole  of  the  ablative  data  (2) 
for  impacts  in  th  excess  energy  regime  and  angles  over  30'  was 


I) 


n 


(16) 


As  more  data  are  accumulated  for  vulcanized  rubber  bonds  between  astrolite  and 
back-up,  much  of  the  scatter  is  being  eliminated.  With  these  data,  the  effects  arising 
from  change  of  backing  material,  incidence  angle,  and  ablative  material  are  more  appar¬ 
ent.  From  the  more  recent  astrolite  data,  the  best  fit  for  hole  sizes  resulting  from  nor¬ 
mal  impacts  is  given  by  Equation  (17)  from  which  the  value  of  ■  for  decomposed  astrolite 
can  also  be  determined 


Substituting  this  value  oi  .  in  Equation  (15)  gives 

•>»  K.,  T„)  ’ 70  (18) 

a  relationship  which  was  derived  for  determining  hole  sizes  only  in  the  low  energy  tran¬ 
sition  regime. 

The  combination  shock  and  ablative-gas  pressure  type  of  ablative  removal  from  an 
approximately  circular  region  following  a  90"  impact  is  only  one  of  several  removal  ac¬ 
tions  observed,  but  is  perhaps  the  most  common  in  the  excess  energy  regime.  This 
regime  is  of  primary  interest  in  the  study  because  it  is  only  in  this  regime  that  signifi¬ 
cant  Impact  damage  can  be  accomplished  behind  the  structure.  Ablative  hole  size  pre¬ 
diction  is  also  very  important  in  the  evaluation  of  thermal  damage  to  components  behind 
the  structure.  This  is  a  complex  task  due  to  the  many  different  ablattves  with  different 
cratering  characteristics  and  types  of  holes  created,  as  shown  in  Figure  15. 

A  case  of  a  low  energy  impact,  insufficient  to  penetrate  the  back-up  is  shown  at  ihe 
top  of  Figure  15.  Noting  the  right  side  of  the  ablative  hole,  it  is  seen  that  holes  were 
opening  up  in  the  laminate  structure  a  little  behind  the  projecting  tips  and  indicate  a  near 
rupture  at  th's  point.  Only  a  slight  increase  in  projectile  energy  would  have  made  con¬ 
siderable  diherence  in  observed  hole  size.  The  center  inset  also  shows  a  rupture  section 
at  the  right  and  underneath  a  thin,  easily  broken  ledge  indicating  possible  onset  of  fracture 
at  this  point.  Both  of  these  cases  show  why  scatter  on  observed  values  of  hole  size  for 
identical  shot  parameters  can  be  plus  or  minus  lO'.Vi.  A  type  of  back-up  failure,  charac¬ 
teristic  of  magnesium  back  ups,  is  shown  in  the  lower  inset.  Large  front  spalls,  shown 
in  the  lower  two  insets,  occur  in  the  low  energy  regime  for  certain  ablative-projectile 
combinations. 

There  are  three  main  energy  regimes  of  ablative  hole  damage  and  are  shown  sche¬ 
matically  in  Figure  16.  The  first  two  sketches  (a  and  b)  show  nonpenetrated  ablative 
structures  with  holes  formed  primarily  by  cratering.  In  the  minimum  energy  case,  the 
hole  is  formed  by  the  intersection  of  the  crater  and  rear  spall,  which  is  usually  pulverized 
and  expelled  by  back-up  rebound.  Sketch  b  shows  the  cratering-spall  hole  created  just 
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Figure  15  -  Photograph  of  the  various  effects  of 
hypervelocity  impacts  on  ablative  structures 
for  complete  performation  of  the  shin  live 
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Figure  16  -  Sketch  of  various  energy  regimes  in  hype rvelocity  impact:  (a)  start 
of  cratering  regime  at  v„  ,  ( b)  beginning  of  transition  regime  at  V,  ,  (c)  model  of 
backup  motion  (luring  transition,  (d)  transition  regime  hole,  (e)  start  of  excess 
energy  regime,  (f)  excess  energy  regime 


before  the  cratering  region  touches  the  back-up.  In  both  cases,  the  high  amplitude  por¬ 
tion  of  the  incident  shock  has  been  relieved  by  side  rarefactions  from  the  projectile  be¬ 
fore  the  shock  can  enter  the  back-up. 

The  two  center  sketches  (c  and  d)  Illustrate  essential  features  of  a  transition  zone, 
where  high  pressure  gas  Is  created  by  a  shock  unattenuated  by  projectile  rarefactions. 
The  decomposition  products  are  confined  to  a  radial  expanding  zone  between  the  ablative 
and  a  slowly  receding  back-up  which  does  not  rupture  or  open  up  until  after  the  gas  pro¬ 
pelled  ablative  fractures  over  a  high  force  gradient  zone.  The  back-up  does  not  acquire 
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high  velocity  because  the  unattenuated  pressure  pulse  traversing  the  back-up  plate  is 
overtaken  in  the  back-up  by  the  side  and  end  rarefactions  from  the  projectile.  This 
regime  is  characterized  uy  very  prominent  bulges  and  cracking  up  to  a  minimum  type  of 
petaling  where  bt,  Dn.  The  resulting  radial  type  of  gas  expansion  has  been  discussed  in 
the  preceding  section  on  theory  and  leads  to  an  energy-thickness  ratio  exponent  of  about 
0,7  in  the  ablative  hole  size  relationship. 

The  sketches  e  and  f  show  the  excess  energy  regime  where  D,,  >  Dn  and  the  back-up 
velocity  was  sufficiently  high  to  insure  rupture  and  approximately  hemispherical  gas  ex¬ 
pansion  before  full,  effective  radial  expansion  of  the  ablative  gases  could  occur.  This 
brief  analysis  then  leads  to  consideration  of  ablative  hole-size  versus  energy-thickness 
ratio  for  the  various  energy  regimes. 

The  ablative  hole  size  versus  energy -thickness  ratio  for  90°  impacts  into  astrolite 
with  steel  back-up  is  shown  in  Figure  17.  First  to  be  considered  is  the  establishment  of 
the  maximum  hole  size  limit  curve  for  low  and  minimum  energy  impacts.  From  Refer¬ 
ence  2,  a  minimum  perforation  energy  was  obtained  for  a  steel  projectile  Impacting 
astrolite  with  a  steel  back-up,  for  ablative  thicknesses  of  1.27  and  2.54  cm.  Under  mini¬ 
mum  perforation  conditions  several  hole  sizes  in  the  ablative  were  examined  for  the  two 
thickness,  This  approximately  established  the  end  of  the  low  energy  or  cratering  regime 
and  the  beginning  of  the  transition  regime,  shown  at  the  upper  right  of  the  shaded  zone  of 
scatter  in  Figure  17.  The  left-hand  limit  point  was  established  by  considering  that  spall 


Figure  17  -  Hole  size  vs  energy-thickness  ratio  for 
steel-backed  astrolite  (90°  impacts) 
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formation  would  complete?  the  ablative  hole  with  impact  energy  sufficient  to  form  a  crater 
depth  about  0.7  of  the  ablative  thickness,  and  hole  size  was  assumed  to  vary  as  <F,  T„)  •  1 
The  high  energy  curve  was  then  drawn  in  to  give  the  best  fit  with  data  for  fully  developed 
petals.  The  transition  curve  was  then  constructed  to  connect  each  of  the  upper  left-hand 
points  (beginning  of  transition  regime)  with  the  experimental  excess  energy  curve.  The 
transition  curve  was  constructed  with  a  slope  determined  by  the  energy  thickness  ratio 
to  the  0.7  power  per  Equation  (18).  It  is  to  be  noted  that  if  the  transition  curves  were 
inappropriately  extended  above  the  excess  energy  regime  intercept,  considerable  error 
would  be  made  in  hole  size  predictions.  The  excess  energy  curve  for  90"  impacts  on 
astrolite  with  aluminum  and  steel  projectiles  is  given  by  Equation  (17)  in  which  hole  size 
is  proportional  to  energy-thickness  ratio  to  the  0.41  power.  There  appears  to  be  one 
curve  within  present  experimental  error  for  ablative  thickness  ratios  of  T„  br,  above  1.5 
in  the  excess  energy  regime  as  contrasted  with  a  curve  for  each  thickness  in  the  low 
energy  and  transition  regimes.  The  formula  for  the  excess  energy  region  for  astrolite 
impacts  at  90  for  energy-thickness  ratios  ranging  from  25  to  125  kj.  cm  gives  hole  size 
predictions  in  good  agreement  with  the  general  ablative  gross  effect  formula  presented 
at  the  Sixth  Symposium  (2). 


Large  Mass  Impacts 

A  four- foot  square  roller-mounted  ablative  target  is  shown  in  Figure  18.  The  abla¬ 
tive  was  4.2-cm  thick  with  a  0.63-cm  thick  steel  rubber-bonded  back-up  mounted  45  to 
the  trajectory.  The  impact  produced  a  35-cm  ablative  hole  and  the  back-up  received  ex¬ 
tensive  petaling  damage  as  shown  by  the  photo  at  the  left.  The  center  edges  of  the  front 
and  rear  photo  show  matching  edges  of  the  target  and  petaling  can  be  seen  to  be  asym¬ 
metric  in  the  direction  of  the  projected  component  of  projectile  motion  along  the  plane  of 
the  back-up. 


ASTROLITE  COMPOSITE 


back  front 

Figure-  1H  -  Massive  impact  (4,1  megaj miles) 
ititu  steel-backed  astrolite  target  at  -IS" 


In  Figure  19  is  presented  a  graphical  summary  of  45'  impacts  into  astrolite  bonded 
to  a  steel  back-up  with  vulcanized  rubber.  The  upper  points  are  the  4.1  and  1.2  megajoule 
impacts  with  respective  366-gm  aluminum  and  100-gm  steel  spheres.  The  lower  points 
each  represent  shots  where  target  and  projectile  dimensions  were  scaled  by  a  factor  of 
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Figure  19  -  Hole  size  vs  energy-thickness  ratio  for 
astrolite  targets  (45*  impact) 


0.16  of  the  corresponding  massive  impacts.  Considering  the  line  for  the  steel  projectiles 
first,  it  is  noted  that  the  intermediate  energy  shots  fell  very  near  the  connecting  line  even 
though  they  represented  impacts  at  a  higher  velocity  (7.8  km/Bec)  and  an  ablative-thickness 
ratio  of  about  2.5.  The  line  connecting  the  massive  and  the  scaled  aluminum  shots  did  not 
fall  upon  the  line  for  the  steel  shots  even  though  they  all  are  in  the  excess  energy  regime. 
This  could  be  attributed  to  a  thin  plate  effect  for  the  low  ablative  thickness  ratio  (Ta  D0  = 
0.68),  or  the  effect  of  change  in  projectile  material  or  a  combination  of  these  effects.  Of 
major  interest  is  the  fact  that  the  hole  diameters  for  these  angled  shots  appear  to  be  pro¬ 
portional  to  the  energy  thickness  ratio  taken  to  about  the  0.54  power.  This  represents  a 
departure  from  the  excess  energy  curve  for  90°  impacts.  One  change  which  would  occur 
in  angled  impacts  would  be  that  multiple  wave  reflections  between  the  back-up  and  the 
metal  projectile  would  no  longer  be  normal  to  each  surface.  Another  effect  is  that  the 
point  of  the  projected  collision  between  the  projectile  and  the  back-up  is  not  at  a  point  of 
maximum  back-up  velocity  from  the  preceding  wave  action. 


Pressure  Instrumentation  in  Ablative  Material 


Though  the  theoretical  model  discussed  previously  is  for  TWPN  impact  by  a  TWPN 
cylinder  at  90°,  the  model  and  methods  outlined  in  the  graphical  wave  representation  can 
be  approximated  to  apply  to  an  impact  such  as  the  large  aluminum  projectile  impacting 
at  45°  and  4.57  km/sec  into  the  astrolite  target.  In  this  experiment,  pressure  data  were 
obtained  using  carbon-resistor  pressure  gages  (10)  placed  in  the  center  of  the  ablative 
about  2,11  cm  in  from  the  ablative  face  and  about  25.4  cm  from  the  impact  center  at  90' 
from  the  trajectory  (see  front  of  target  in  Figure  18).  To  properly  evaluate  the  pressure 
recorded,  reference  is  made  to  the  model  shown  in  Figure  8.  Had  the  projectile  in  this 
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figure  been  of  a  higher  specific  acoustical  impedance,  an  increase  of  pressure  would 
have  occurred  when  the  reflected  shock  wave  component  from  the  back-up  interacted 
with  the  projectile  material  interface.  This  reflected  component  then  would  have  been 
propagated  in  a  spherical  manner  back  toward  the  steel  back-up.  When  this  projectile 
reflected  shock  component  reached  the  steel  back-up,  a  second  shock  component  reflec¬ 
tion  will  occur  and  will  propagate  back  into  the  ablative,  The  actual  massive  impact 
problem  is  further  complicated  by  the  45"  angle  of  attack  so  that  no  attempt  is  made  to 
predict  pressure  wave  values  and  the  effect  of  all  possible  ray  paths  of  propagation.  The 
pressure-time  record  of  these  waves  in  the  massive  impact  experiment  does  show  pres¬ 
sure  magnitude  and  arrival  time  of  these  components.  In  Figure  20,  P,  represents  the 
maximum  value  of  the  direct  pressure  wave  from  the  impact  already  under  the  influence 
of  the  rarefaction  ray  from  the  surface.  Rarefaction  attenuation  continued  to  predominate 
until  shock  reflection  from  the  back-up  started  to  act  at  P2  just  0.405  ,  sec  after  the  main 
wave  arrival  or  scope  triggering.  At  about  1.1  .  see  after  the  scope  triggered,  signals 
from  the  projectile  reflection  of  the  back-up  reflection  wave  started  to  act  and  reduce  the 
rarefaction  effects.  Pressure  rose  slightly  to  a  maximum  value  P4  before  surface  rare¬ 
factions  again  predominated.  The  final  compression  wave  from  the  back-up  reflection  of 
the  projectile  reflection  component  reached  a  maximum  value  of  p  . 


CONCLUSIONS 

It  has  been  shown  that  a  simple  two-dimensional  wave  analysis  can  be  applied  to  gain 
a  good  idea  of  momentum  transfer  from  the  projectile  to  the  target.  Augmentation  ratios 
predicted  are  consistent  with  observed  momentum  augmentation,  and  radial  momentum 
determined  from  behind  the  target  cratering  studies  (1,9)  are  in  line  with  early  stage 
predictions  based  on  wave  analysis.  Future  work  is  necessary  to  adapt  wave  analysis 
techniques  to  impacts  made  at  incidence  angles  less  than  90°,  and  experimental  pressure¬ 
time  recordings  can  do  much  to  guide  the  theoretical  approach  into  these  more  realistic 
modes  of  attack. 


Figure  .JO  -  Pressure-time  curve  for  a 
massive  astrolite  impact  (.15.4  cm  from 
impact  point) 
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Three  regimes  of  ablative  hole  damage  exist:  low  energy,  transition,  and  excess 
energy.  In  the  transition  "egime  there  is  a  separate  curve  for  each  ablative  thickness, 
and  the  ablative  hole  size  varies  as  the  energy-thickness  ratio  to  the  0.7  power.  Inap¬ 
propriate  extension  of  this  curve  beyond  the  point  where  high  energy  transition  occurs 
can  lead  to  large  errors  in  hole  size  prediction.  The  excess  energy  regime  is  of  greatest 
interest  because  of  the  impact  damage  potential  to  components  located  behind  the  target. 
The  curve  for  90°  impacts  into  astrolite  with  steel  back-up  show  that  the  hole  size  is 
proportional  to  the  energy-thickness  ratio  to  the  0.41  power  in  the  thickness  range  where 
2.8  s  Tn  D0  >  1.4.  Values  of  hole  size  in  the  25  to  125  kilojoules/cm  range  are  very  close 
to  those  presented  in  Reference  2  for  impacts  into  ablatives  in  the  excess  energy  regime. 

The  experimental  shots  made  at  45°  have  shown  that  hole  size  is  proportional  to 
energy-thickness  ratio  to  the  0.54  power.  This  increase  in  power  could  be  due  to  a 
change  from  the  normal  wave  interaction  between  back-up  and  projectile  which  occurs 
in  90°  impacts.  Hole  sizes  from  shots  made  with  aluminum  projectiles  at  Ta/D0  <”  0.7 
fell  below  the  curve  for  ratios  greater  than  1.4,  indicating  a  thin  plate  effect. 
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diameter  of  hole  in  ablative 

diameter  of  hole  in  back-up 

diameter  of  projectile 

energy  in  impacting  projectile 

force 

acceleration  of  gravity 

constant  dependent  upon  target  and  projectile  material,  X  =  3y-2 
constants  empirically  determined 

constant  empirically  determined  (dependent  upon  target  and 
projectile  material)  a  -  l/3y-l 

constant  empirically  determined  (low  energy  transition) 

H  =  \/2y  -  1 


M0  gm-km/sec 
M;  gm-km/sec 
M.  r  gm-km/sec 

n 


P  KB 
Pc  KB 
PB  KB 
AP  j  KB 
t  n  sec 
U  km/sec 


momentum  of  impacting  projectiles 
axial  momentum  of  ith  ring  of  particles 
radial  momentum  of  ith  ring  of  particles 
number  of  waves  passed  at  chosen  time  and  place 
pressure  in  material  at  given  location  and  time 
initial  pressure  in  material 
pressure  of  gas  in  ablative  cavity 
incremental  change  in  pressure  due  to  ith  wave 
time 

shock  wave  velocity  (reference  to  undisturbed  region) 
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u  km/ see  particle  velocity  in  material  at  given  location  and  time 

u  km/sec  initial  velocity  of  material  or  projectile 


11 .. 

km/sec 

velocity  of  projectile  and  target  ejecta 

11  L 

km/ sec 

axial  particle  velocity  of  ith  region 

11  r 

km/sec 

radial  particle  velocity  of  i  th  region 

'u,l 

km/ sec 

. ? 

incremental  change  in  particle  velocity  due  to  ith  wave 

....1 . r  . . . 

’ll 

W 

gm/ sec 

rate  of  material  ejection 

*1 

gm 

weight  in  i  th  wave  region 

ratio  of  specific  heat  (decomposed  ablative  gas' 
ratio  of  densities,  .  .  0 


incremental  change  in  .. 

gm/cm3  density  of  shocked  material  under  pressure 

distance  projectile  advances  into  target  divided  by  D„  diam¬ 
eter  (assuming  u0  constant) 

cm  thickness  (depends  on  scale  factors  of  the  experiment) 
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SUMMARY 

Work  has  continued  under  the  Aerothermal  Phase  of  the  program  with  em¬ 
phasis  on  the  thermal  performance  of  foam -filled  axisymmetric  models  having 
simulated  hypervelocity  impact  perforations. 

The  Malta  Rocket  Exhaust  Internal  Heudng  Test  Program  was  concluded 
with  the  completion  of  Tests  Number  6  and  7.  The  principal  achievements  of  this 
program  were  the  establishment  of  a  correlation  for  the  rate  of  energy  influx 
through  a  heat  shield  perforation  and  the  determination  of  the  thermal  performance 
of  rigid  urethane  foam  within  perforated  models,  detailed  results  of  Tests  Number 
fi  and  7  are  presented  herein,  together  with  a  summary  of  the  overall  program  and 
principal  results.  The  application  of  the  correlation  to  a  high  performance  ICBM 
re-entry  vehicle  indicates  that  an  internal  heating  load  of  about  200  btu  per  square 
inch  of  perforation  area  is  accumulated  by  50, 000  feet  altitude  for  a  perforation 
occurring  at  100, 000  feet  altitude.  The  particular  foams  tested  showed  an  energy* 
absorbing  capability  of  about  2000  and  3000  btu/lb  for  densities  of  3.0  and  7.8  lb/ft  , 
respectively. 

Two  new  tasks  were  added  to  the  program  in  August.  Under  the  first  task, 
an  extensive  investigation  of  the  thermal  performance  of  lightweight  foams  having 
simulated  hypervelocity  impact  damage  will  be  conducted.  Axisymmetric  foamed- 
filled  models  will  be  tested  in  the  GE-RSD  5  -  Megawatt  Air  Arc  and  Malta  Rocket 
Exhaust  facilities.  The  choice  of  foams  will  be  based  on  screening  tests  performed 
in  a  blowtorch  stream.  Calorimeter  models  have  been  fabricated,  and  the  arc 
shroud  model  configuration  is  being  designed.  A  tentative  list  of  promising  foam 
candidates  has  been  formulated. 

Under  the  second  new  task,  a  thermostructural  kill  study  is  being  conducted, 
with  the  objectives  of  (1)  formulating  a  simplified  computational  model  for  thermo¬ 
structural  kill  of  ICBM  re-entry  vehicles  due  to  hypervelocity  impact  perforations 
and  (2)  applying  the  computational  model  to  selected  advanced  re-entry  vehicle  de¬ 
signs  of  two  general  classes:  (a)  hardened  and  (l))  unhardened.  Initial  study  under 
this  task  has  resulted  in  the  selection  of  seven  vehicle  configurations  having  bal¬ 
listic  coefficients  from  1700  to  4400  Ib/ft^. 
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I.  MALTA  ROCKET  EXHAUST  INTERNAL  HEATING  TEST  PROGRAM 
A.  General 


The  Malta  Rocket  Exhaust  Internal  Heating  Test  Program  was  completed 
with  the  performance  of  Tests  Number  6  and  7.  The  general  objectives  of  the 
program  and  detailed  results  of  Tests  Number  1  through  5  were  presented  in  the 
preceding  semiannual  progress  report  (Ref.  1).  A  condensed  summary  report  of 
the  entire  program  was  also  published  (Ref.  2). 

The  test  agenda  and  objectives  of  the  program  are  summarized  in  Table  1. 
The  tests  were  conducted  on  identical  sphere-cone  models  of  9.8  inches  base  diam¬ 
eter  consisting  of  a  laminated  phenolic  nylon  heat  sliicld  hard-lwndcd  to  an  Inconel 
structure.  A  single  perforation  was  located  on  the  conical  skirt  of  each  model. 
Tests  Number  1  through  4  featured  open  internal  compartments,  while  Tests  Num¬ 
ber  5  through  7  feat  tred  foam  -filled  internal  compartments.  In  Tests  Nurn!>er  1 
through  5,  30  thermocouples  were  used  to  measure  heat  flux  distribution  to  the 
structure.  For  Tests  6  and  7,  instrumentation  was  modified  to  include  10  thermo¬ 
couple  probes  positioned  within  the  foam  to  indicate  the  foam  surface  recession 
history. 

TABLE  1.  MALTA  ROCKET  EXHAUST  INTERNAL  HEATING  TEST  SUMMARY 
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Nominal  test  conditions  for  all  tests  are  shown  in  Table  2. 
TABLE  2.  MALTA  PIT  FOUK  TEST  CONDITIONS 


Free-stroam  Mach  num tier 

3.0 

Model  stagnation  conditions: 

Pressure 

170  psia 

Enthalpy 

3250  btu/lb 

Temperature 

5600  °H 

Cone  static  pressure 

22  psia 

Test  medium 

Exhaust  gases 
(y  -  1.2,  mol,  wt.  -  25) 

The  model  for  Test  Number  7  is  shown  mounted  on  the  pit  number  four  test  stand 
in  Figure  1,  prior  to  gimbaling  the  rocket  exhaust  nozzle  into  position  for  test. 


Figure  1.  Model  on  Test  Stand  (Malta  Test  No.  7) 
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B.  Results  of  Tests  Number  6  imd  7 


1.  General 


The  principal  objective  of  Tests  No.  G  and  7  was  to  determine  the  heat 
protection  characteristics  of  a  Lightweight  urethane  foam  within  axisymmetrlc 
models  containing  simulated  hypervelocity  impact  damage.  In  each  test,  the  foam 
used  was  a  3  lb/ft^  rigid  urethane  foam  (Hetrofoam  368)  made  by  the  Hooker 
Chemical  Co. 

To  accomplish  this  objective,  thermocouple  beads  were  mounted  at  the 
end  of  0,  2!>- inch  diameter  probes  m.ulo  of  GLv  Series  100  ablation  material  which 
were  installed  within  the  foam  at  ten  locations  as  shown  in  Figure's  2  and  3.  The 
probes  were  threaded  into  the  Inconel  structure  and  were  arranged  to  minimize 
interference  between  adjacent  probes. 

Cross-sections  of  the  simulated  impact  cavity  for  each  test  are  also 
shown  in  Figures  2  and  3.  These  cavities  were  idealized  representations  of 
hyperveloeity  impact  damage  into  phenolic  nvlon-steel  targets  backed  by  urethane 
foam  (Hetrofoam  3G8),  The  impart  tests  were  performed  by  the  Naval  Research 
Laboratory. 

The  tests  were  intentionally  run  for  a  long  enough  time  (15  seconds)  to 
completely  remove  the  foam  during  the  test,  thereby  causing  heating  of  the  Inconel 
structure  following  foam  removal  by  ablation.  Structure  temperatures  were 
monitored  by  20  thermocouples  located  on  the  structure  at  the  shield-structure 
interlace,  Four  static  pressure  taps  were  used  to  measure  pressure  within 
each  model.  Thermocouple  and  pressure  tap  locations  are  shown  in  Figures  2 
anti  3, 

Pre-test  and  post-test  photographs  of  the  models  for  Tests  No.  6  and  7  are 
given  in  Figures  4  through  7. 

Figure  7  shows  a  longitudinal  crack  0,  25  inch  in  width  that  developed 
downstream  of  the  perforation  in  Test  No.  7.  Pastes  movies  showed  that  the 
crack  occurred  at  12  seconds  during  the  run,  A  similar  crack  developed  during 
Test  No.  2  (see  Ref,  1).  In  Figure  5,  scale  dip.  sits  can  be  seen  on  the  backplate 
within  the  model.  Similar  deposits  were  present  in  Model  No.  7. 

Typical  probe  thermocouple  temperature  histories  for  Tests  No.  6  and  7  are 
shown  in  Figures  8  and  9,  Due  to  the  low  thermal  conductivity  of  the  foam,  no 
significant  temperature  response  occurred  until  the  foam  surface  receded  to  the 
thermocouple  location.  The  thermoc. .titles  thus  functioned  a«  ablution  sensors 
by  locating  the  foam  surface  at  various  times  during  the  test.  After  the  foam 
surface  receded  past  the  thermocouples,  their  readings  appea*  to  level  off  in 
approximate  equilibrium  with  the  local  gas  recovery  temperature  at  1500  to 
2200°  F. 


Plots  of  foam  surface  recession  vs.  time  as  determined  from  probe 
thermocouple  records  are  given  in  Figures  10  and  11  for  Tests  No,  6  and  7. 
In  several  cases,  the  final  surface  location  for  a  given  thermocouple  ray  was 
obtained  from  the  response  of  thermocouples  mounted  on  the  structure  in  line 
with  the  probe  thermocouples. 
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Figure  2.  Cross-Section  of  Test  Model  No. 


Figure  3.  Cross-Section  of  Test  Model  No. 
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Figure  4.  Malta  Test  No.  6  -  Pre-Test  View 
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Figure  6.  Malta  Test  No.  7  -  Pre-Test  View 
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Figure  7.  Malta  Test  No.  7  IHsl -Test  View 
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Figure  8,  Typical  Temperature  Histories  -  Malta  Test  No.  6,  Ray  B 
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2.  Thermal  Performance  of  Hetrofoam  368 

Figures  10  and  11  were  used  to  construct  approximate  foam  surface 
recession  contours  as  shown  in  Figure  12  and  13.  These  contours  were  used 
to  approximate  the  weight  of  foam  removed  ul  various  times  during  the  test.  The 
foam  cavity  cross-section  in  a  plane  perpendicular  to  the  cavity  axis  was  assumed 
to  be  elliptical,  with  a  variation  of  b/a  (I)  major  axis,  a  -  minor  axis)  corres¬ 
ponding  to  the  final  foam  cavity  shape  of  Test  No.  3, 

The  thermal  performance  of  the  foam  was  characterized  by  an  effective 
heat  of  ablation  Q*  defined  as: 


in  which  aE  total  energy  entering  perforation  cavity  and 
AM  ■  weight  of  foam  removed  by  ablation. 

Values  of  AM  at  various  times  during  the  run  were  determined  as  described  in 
the  preceding  paragraph.  Values  of  aE  Were  determined  by  assuming  a  constant 
rate  of  energy  addition  of  133  btu  sec  to  the  cavity,  based  on  the  result  obtained 
in  Test  No.  2  for  the  same  size  perforation.  Therefore  equation  (1)  is  an  arbitrary 
definition  for  Q*  in  which  the  effects  of  perforation  area  enlargement  and  foam 
decomposition  gases  generated  during  the  test  are  included  in  the  definition. 

The  results  of  the  analysis  described  above  are  given  in  Table  3.  The 
result  given  for  Q*  for  Test  No.  5  supersedes  the  value  given  in  Ref.  1,  which 
was  incorrect  due  to  an  error  in  the  calculation  of  ablated  foam  volume.  Con¬ 
sistent  results  were  obtained  from  computations  performed  for  times  of  2  and  4 
seconds,  using  the  contours  of  Figures  12  and  13. 


TABLE  3.  VALUES  OF  FOAM  HEAT  OF  ABLATION  Q*  FOR 
HETROFOAM  368 


TEST  NO. 

5 

6 

7 

Foam  density  (lb/ft®) 

7.  8 

3.0 

3.0 

Q*  (btu/lb) 

2900 

2300 

1800 

These  values  of  Q*  are  subject  to  uncertainties  on  the  order  of  ±25  percent  due 
to  uncertainties  in  the  exact  values  of  aE  and  AM. 

In  Tests  No.  6  and  7,  the  structure  heat  flux  distribution  following  ablation 
of  the  foam  indicated  that  the  maximum  heating  occurred  downstream  of  the  per¬ 
foration,  rather  than  opposite  the  perforation  as  in  Tests  1  through  4.  This  shift 
in  internal  heat  flux  distribution  is  attributed  to  the  thinner  downstream  edge 
caused  by  ablation  during  the  longer  runs.  Exact  heat  fluxes  and  rate  of  energy 
absorption  could  not  be  determined  following  foam  ablation  due  to  the  deposition 
of  scale  formed  from  the  melted  Inconel  structure  adjacent  to  the  perforation. 
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(Sot*  Figure  5.)  The  scale  had  a  relatively  thick,  c  rusty  appearance,  and  con¬ 
stituted  a  significant  added  thermal  resistance  which  reduced  the  local  heat  flux 
due  to  increased  surface  temperatures.  This  effect  was  verified  experimentally 
by  subjecting  scaled  and  unsealed  areas  of  the  backplatc  to  identical  heating  con¬ 
ditions  from  a  blowtorch.  The  temperature  response  of  a  thermocouple  on  the 
back-face  of  the  scaled  area  Indicated  a  heat  flux  of  approximately  half  that  for  the 
unsealed  area. 

3.  Perforation  Enlargement  Due  to  Ablation 


Tracings  and  measurements  of  the  perforation  contours  were  made 
directly  from  the  models  after  testing  for  Models  No.  4  through  7.  Pro-test 
and  post-test  contoui  s  of  the  perforations  in  the  heat  shield  and  structure  are 
shown  in  Figure  14.  More  heat  shield  material  was  removed  than  metal  structure, 
despite  the  greater  heat  of  ablution  of  phenolic  nylon  relative  to  that  of  Inconel. 
This  effect  is  due  partly  to  the  higher  heat  flux  near  the  upper  corner,  but  may 
also  be  the  result  of  poorer  ablation  performance  due  to  high  vorticlty  in  the  re¬ 
attachment  region  of  the  separated  shear  layer  which  develops  over  the  perfora¬ 
tion.  The  steep  gradient  in  ablation  of  the  downstream  edge  of  the  perforation 
is  dtown  in  the  sketch,  which  is  typical  of  results  for  Tests  No.  5  through  7, 
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It  is  seen  from  Figure  14  that  significant  perforation  enlargement  due  to 
ablation  occurred  for  the  longer  test  durations  of  Tests  No.  6  and  7.  Computa¬ 
tions  for  these  cases  yielded  the  results  shown  in  Table  4. 

TABLE  4.  PERFORATION  ENLARGEMENT  DUE  TO  ABLATION 


TEST  NO. 

6 

7 

Percent  Increase  in  heat 
shield  perforation  area 

43 

37 

Percent  increase  in  Inconel 
liner  perforation  area 

19 

18 

4.  Pressure  Measurements^ 

The  pressure  measurements  for  Test  No,  7  are  plotted  in  Figure  15. 
Little  pressure  rise  was  observed  for  tap  locations  on  the  backplate,  due  to 
clogging  by  pieces  of  foam.  The  taps  located  directly  beneath  the  upstream 
edge  of  the  perforation  recorded  a  pressure  history  quite  similar  to  that  which 
was  obtained  for  Test  No.  2,  Thus,  it  appears  that  the  presence  of  foam  de¬ 
composition  gases  within  the  model  has  no  significant  effect  on  the  pressure 
history. 
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RUN  TIME  ~  SECONDS 

Figure  15.  Internal  Pressure  Histories  -  Malta  Test  No.  7 
5.  Foam  Characteristic s 

The  results  of  Teats  No.  5,  6,  and  7  shown  in  Table  3  indicated  relatively 
good  performance  of  Hetrofoam  368  for  the  Malta  rocket  exhaust  environment. 
This  observation  prompted  additional  laboratory  investigations  of  the  structural 
and  thermal  degradation  characteristics  of  this  foam. 

Compressive  tests  of  1  inch  cubes  were  performed  on  an  Instron  universal 
testing  machine  to  determine  the  modulus  of  elasticity  and  ultimate  stress  in  the 
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foam-rise  direction  and  perpendicular  to  it.  The  test  results  are  shown  in  Table  5, 
and  indicate  higher  strength  in  the  foaming  direction,  as  well  as  increased  strength 
with  increased  density. 


TABLE  5.  COMPRESSIVE  TEST  RESULTS  FOR  HETROFOAM  368 


Density 

(lb/ft3) 

Modulus  of 

Ultimate 

Elasticity 

(psi) 

Strength 

(psi) 

Loading  Condition 

3.0 

1800 

60 

Direction  of  foaming 

3.0 

700 

52 

Perpendicular  lu  direction  of  foaming 

7.3 

4900 

152 

Direction  of  foaming 

7.3 

_ 

3300 

115 

Perpendicular  to  direction  of  foaming 

Thermal  degradation  tests  were  performed  by  heating  rectangular  samples 
to  1000°F  by  means  of  a  Nichrome  wire  connected  to  a  Variac.  In  the  first  test, 
the  temperature  was  increased  slowly  (from  room  temperature  to  1Q00°F  took 
2  hours);  in  the  second  test,  a  more  rapid  heating  rate  was  used  to  reach  1000°F 
in  a  few  minutes.  Samples  of  both  densities  (3. 0  and  7. 3  lb/ft3)  reacted  similarly; 
they  warped,  discolored,  charred,  and  finally  broke  due  to  spring  tension  induced 
by  the  heater  wire.  Charring  began  at  a  temperature  of  approximately  500°F. 

The  foam  smoked  considerably  in  the  higher  heating  rate  test,  but  no  sign  of  com¬ 
bustion  was  observed.  It  was  noted  that  the  foam  burned  in  the  flame  of  a  match, 
but  was  self -extinguishing  upon  the  formation  of  a  char  layer  on  its  surface. 


C.  Final  Correlation  of  Internal  Heating  to  Open  Volumes 

In  Ref.  1,  a  semi-empirical  correlation  of  wind  tunnel  and  rocket  exhaust 
test  results  was  developed  for  the  rate  of  energy  absorption  by  the  walls  of  axi- 
symmetric  models  having  single  perforations.  The  correlating  parameter  was 
based  on  the  compressible  turbulent  shear  layer  analysis  of  Chow  and  Korst  (Ref.  3), 
and  included  the  Crocco  number  C  and  a  mixing  similarity  parameter  o. 

Further  study  of  the  results  showed  that  the  effect  of  C  and  a  can  be  in¬ 
corporated  into  the  empirical  proportionality  constant  with  little  loss  in  accuracy, 
due  to  the  compensating  effect  of  Mach  number  on  these  terms.  The  resulting  cor¬ 
relation  is  shown  in  Figure  16,  and  is  considerably  simpler  to  evaluate  than  the 
previous  correlation.  The  correlating  parameter  used  in  Figure  16  is  equivalent 
to  the  simple  result  obtained  by  Donaldson  (Ref.  4)  by  means  of  a  turbulent  shear 
layer  mixing  analysis,  with  the  constant  adjusted  to  agree  with  the  experimental 
results.  The  success  of  this  method  of  correlation  over  a  range  of  3  decades  in 
the  correlating  parameter  and  for  two  widely  different  test  environments  lends 
credence  to  the  shear  layer  model  for  the  case  of  a  single  perforation. 

The  resultant  correlation  equation  for  the  rate  of  energy  absorption  by  the  struc¬ 
ture  is: 

dE/ dfl  =  ■  0105  y  ^  —  (2) 
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with  symbols  being  defined  in  Figure  16.  The  scatter  of  the  bevelled  and  impact 
perforation  results  about  the  mean  correlation  line  for  cylindrical  perforations 
indicates  that  the  correlation  is  a  reasonable  representation  for  actual  impact 
perforations. 


D.  Application  of  Results  to  ICDM  Re-entry 


Equation  (2)  was  applied  to  the  C-l  Advanced  Target  Vehicle  design  of  Ref.  5. 
This  hypothetical  re-entry  vehicle  is  a  slightly  blunted  cone  with  the  following 
characteristics: 


RntHstic  coefficient: 
Semi-vertex  angle: 
Re-entry  velocity: 
Re-entry  path  angle: 


3000  lb/ft2 
11  degrees 
25, 000  feet/second 
23  degrees 


Figure  16,  Correlation  of  Internal  Heating  Results 
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The  results  are  shown  In  Table  6,  and  apply  for  a  single  perforation  anywhere  on 
the  cone,  assuming  intercept  to  occur  at  an  attitude  of  100, 000  feet. 

TABLE  6.  ENERGY  INFLUX  THROUGH  PERFORATION  FOR  C-l  VEHICLE 


Z  -Altitude(ft. ) 

70,000 

50, 000 

30, 000 

E/A  integrated  energy  pm*  unit 
perforation  area,  (btu/inz)  from 
100,000  ft.  to  Z  ft. 

50 

170 

375 

These  results  neglect  the  effect  of  perforation  enlargement  during  re-entry.  From 
calculations  of  this  type,  lethal  perforation  areas  were  presented  in  Ref.  1  for 
various  structural  designs  of  the  C-l  class  of  vehicle  as  a  function  of  vehicle  size, 
intercept  altitude,  and  kill  altitude. 

Estimates  of  weight  penalties  to  a  given  vehicle  design  by  the  use  of  internal 
foam  for  hardening  against  thermal  effects  can  be  obtained  by  use  of  a  foam  heat 
of  ablation  on  the  order  of  2000  to  3000  btu/lb,  as  obtained  in  Tests  5  through  7.  * 
However,  choice  of  foam  density  must  be  related  to  foam  compressive  strength 
considerations.  Cone  pressures  for  high  ballistic  coefficient  vehicles  will  be 
several  times  higher  than  that  of  the  Malta  environment,  such  that  compressive 
collapse  of  the  lighter  foams  may  occur. 


^This  result  Is  strictly  applicable  only  for  conditions  simulated  by  the  Malta  test 
environment.  Higher  values  of  Q*  would  be  expected  for  higher  enthalpy  levels 
of  flight.  This  effect  will  be  pursued  in  a  new  phase  of  the  program,  as  discussed 


in  Section  II  which  follows. 


SECRET 


H-22 


SECRET 


n.  THERMAL  PERFORMANCE  OF  FOAMS 

A  one-year  study  was  begun  in  August  1964  with  the  objective  of  developing 
an  understanding  of  the  heat  protection  characteristics  of  lightweight  foam  having 
simulated  hypervelocity  impact  damage. 

The  use  of  lightweight  foams  has  been  proposed  as  a  means  of  heat  pro¬ 
tection  for  the  interior  of  a  re-entry  vehicle  whose  shield/ structure  may  be 
perforated  by  hypervelocity  impact.  Since  impact  tests  have  demonstrated  that 
foam  fillers  placed  behind  outer  shells  are  effective  in  absorbing  the  impact  energy, 
the  use  of  internal  foam  would  appear  to  have  a  dual  purpose. 

Li  order  l"  avoid  excessive  weight  penalties,  foam  densities  of  less  than 
15  lb/ft3  are  desirable.  Essentially,  no  thermal  performance  data  for  this  class 
of  material  exist  at  nresent. 

The  study  is  being  conducted  under  two  separate  but  related  phases: 

Phase  1.  Study  if  Fundamental  Ablation  Process  for  Cavity  Heating 

Phase  ",  Response  of  Selected  Foam: 

In  Phase  1,  effort  will  be  ■.  oncenlrated  on  the  understanding  of  the  basic 
ablation  phenomena  associated  with  cavity  heating.  In  Phase  2,  the  response  of 
selected  foams  having  simulated  impact  cavities  will  be  determined  in  GE-RSD's 
5-Megawatt  Air  Arc  Facility.  A  detailed  description  of  each  phase  of  the  program 
is  given  below . 

A.  Program  Plan 

1.  Phase  1.  Study  of  Fundamental  Ablation  Process  for  Cavity  Heating 

In  tills  phase  of  the  program,  experimental  studies  will  be  conducted 
according  to  the  test  program  shown  in  Table  7,  using  the  model  configurations 
shown  in  Figure  17.  Emphasis  will  be  placed  on  developing  an  understanding  of  the 
basic  ablation  process  for  lightweight  foams  having  simulated  impact  cavities.  A 
discussion  of  each  group  of  tests  follows: 

Tests  1  and  2:  Calorimeter  models  will  be  tested  to  determine  the 
calorimeter  heat  flux  distribution  within  various  shapes  of  cavities. 

External  pi  ossure  measurements  will  be  made  from  which  local  flow 
conditions  approaching  the  cavity  will  be  determined. 

Tests  3  through  8:  Foam  response  for  two  densities  of  urethane  foam 
will  be  determined  in  three  different  facilities  in  order  to  determine  the 
effects  of  widely  varying  test  environments.  The  characteristics  of  these 
facilities  are  shown  in  Table  8, 

Tests  9  and  10:  The  effect  of  simulated  impact  cavity  size  and  shape  on 
foam  performance  will  be  determined  for  a  given  foam  in  a  given  test 
environment.  The  foam  degradation  will  be  correlated  with  the  calorimeter 
heat  flux  distrilmtion  of  Test  1  for  each  cavity  shape. 
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SHROUD 


(a)  ARC  SHROUD  CONFIGURATION 

(SAME  MODEL  WILL  BE  USED  FOR  MALTA  PIT  ONE  TESTS) 


(b)  ARC  TUNNEL  CONFIGURATION 
Figure  17.  Thermal  Test  Model  Configurations 
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TABLE  7.  THERMAL  TEST  PROGRAM  -  PHASE  ONE 


Test 

No. 

Foam 

Density 

(lb/ft3) 

Impact 

Cavity 

No. 

Test 

Facility 

Test 

Objective 

1 

1, 

2,3 

Arc  Shroud 

Model  Calibration 
(Calorimeter  heat  flux 
pressures). 

2 

1. 

2,3 

Arc  Tunnel 

S 

4 

5 

Ur  el 

luino 

7. 

5 

Arc  Shroud 

Arc  Tunnel 

Malta  Pit  One 

Determine  effect  of 
test  environment  on 
foam  response  for  two 
densities  of  urethane 
foam  having  a  given 
impact  cavity. 

6 

3. 

0 

Arc  Shroud 

7 

Arc  Tunnel 

8 

Malta  Pit  One 

g 

10 

2 

3 

Arc  Shroud 

1 

_ 

Determine  effect  of 
impact  cavity  size  and 
shape  on  foam  response 
(comparison  includes 
Test  No.  6). 

EXAMPLES  OF  IMPACT  CAVITIES 


TABLE  8.  THERMAL  TEST  FACILITY  CHARACTERISTICS 


Facility 

GE-RSD 

5  MW 

Arc  Tunnel 

GE-RSD 

5  MW 

Are  Shroud 

Malta  Rocket 
Pit  One 

Stagnation  enthalpy 
(btu/lb) 

6000  -  20,  000 

6000  -  20,  000 

3250 

Model  Stag,  press,  (atm) 

0.  1  -  ).  0 

1  -  5 

8 

Free  Stream  Mach  No. 

5  -  7 

Subsonic 

2.  5 
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2.  Phase  2.  Response  of  Selected  Foams 

A  review  of  available  foams  will  be  made  to  select  several  possible  materials 
suitable  for  this  application.  Based  on  present  knowledge,  it  is  expected  that  foams 
of  the  following  classes  will  be  considered:  (a)  epoxy,  (b)  silicone,  (c)  urethane, 

(d)  phenolic,  (e)  ceramic,  and  (f)  composites  such  as  ceramic  backed  with  urethane. 
In  addit’on,  consideration  will  be  given  to  lightweight  materials  such  as  balsawood 
(fire-retardant  impregnated)  and  Insulcork. 

The  capacity  of  a  foamed  solid  to  serve  as  a  thermal  barrier  is  a  function 
of  a  number  of  interrelated  variables,  including:  (a)  composition  of  material, 

(b)  wall  size,  (c)  pore  size,  (d)  pore  size  distribution,  and  (e)  flammability  of  the 
foam.  The  materials  chosen  for  thermal  evaluation  will  provide  a  broad  spectrum 
in  anticipated  structure  and  potential  thermal  performance.  As  an  example,  one 
or  more  specific  formulations  will  be  varied  to  produce  a  range  of  densities,  with 
corresponding  variations  in  structural,  mechanical  and  thermal  properties. 

Promising  foam  candidates  will  be  evaluated  in  a  thermal  test  program  as 
outlined  in  Table  9.  An  initial  screening  phase  will  be  conducted  in  which  flat 
samples  will  be  subjected  to  a  blowtorch  test  to  determine  relative  performance. 

The  more  promising  foams  will  be  tested  in  axisymmetric  model  configurations 
(Figure  17)  according  to  the  test  schedule  of  Table  9.  Simulated  impact  cavities 
will  be  provided  in  the  foam.  A  discussion  of  each  group  of  tests  follows: 

Tests  1  -  20:  Flat  circular  disk  specimens  of  foam  will  be  screened 
for  thermal  response  in  a  blowtorch  stream. 

Tests  21  -  26:  The  relative  performance  of  several  foams  will  be 
determined  for  a  given  simulated  impact  cavity  in  a  given  test 
environment. 


TABLE  9.  THERMAL  TEST  PROGRAM  -  PHASE  TWO 


Test 

Foam 

Density 

Impact 

Test 

Test 

No. 

(Choices  are  for 

lb/ft3 

Cavity 

Facility 

Objective 

illustrative 

No. 

purpose  only) 

1  -  20 

Miscellaneous 

3  -  15 

None 

Blowtorch  „ 

Screen  candidate 

q  =  20  btu/ft 

foams  for  further 

sec 

testing.  Flat  test 

P  =  1  atm. 

samples. 

21 

Epoxy 

3 

1 

Arc  Shroud 

Determine  relative 

performance  of 

22 

Epoxy 

7 

several  foams  hav¬ 
ing  a  given  impact 

23 

Silicone 

10 

cavity  ---  (com¬ 
parison  includes 

24 

Ceramic 

10 

Tests  No.  3  and  6 
of  Table  7). 

25 

Composite 

10 

26 

Phenolic 

3 
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In  all  tests,  thermocouples  will  be  placed  at  selected  locations  to  monitor 
the  heat  flux  reaching  the  metal  substructure  due  to  Internal  heating.  Final  foam 
degradation  will  be  measured  and  correlated  with  the  measured  calorimeter  heat 
flux  distribution  for  the  given  cavity  shape  to  determine  the  effective  heat  of  abla¬ 
tion  for  the  foam.  In  addition,  arrays  of  thermocouples  will  be  embedded  within 
the  foam  on  selected  tests  to  obtain  approximate  recession  rate  histories. 

D.  Progress  to  Date 


1.  Phase  1 

The  arc  shroud  calorimeter  model  has  been  designed  and  fabricated  and 
is  shown  in  Figure  18.  This  model  will  lie  used  to  determine  the  convective  heat 
flux  distribution  to  three  simulated  impact  cavities.  In  Figure  18,  cavities  num¬ 
ber  1  and  2  are  geometrically  similar  to  check  scale  effects,  while  cavity  number 
3  represents  the  anticipated  effect  of  an  oblique  impact.  The  cavities  have  been 
located  fairly  far  back  on  the  model  in  order  to  facilitate  transition  to  a  turbulent 
boundary  layer  upstream  of  the  cauiles.  The  heat  flux  distributions  obtained  in 
the  calorimeter  model  test  will  be  used  in  the  analysis  of  the  ablution  performance 
of  foams  having  the  same  initial  impact  cavity  shapes. 

The  basic  arc  shroud  foam  test  models  have  been  designed.  The  external 
configuration  for  the  models  will  be  the  same  as  the  calorimeter  model  of  Figure 
18,  including  a  simulated  Impact  cavity;  however,  the  internal  portion  of  the 
models  will  be  hollow  to  contain  various  foams  for  thermal  evaluation,  The  model 
shell  will  consist  of  a  one-piece  phenolic  nylon  heat  shield  of  0.  25  inch  thickness 
bonded  to  an  0.  080  inch  aluminum  strui  lure.  A  graphite  shroud  will  channel  the 
arc  flow  over  the  model  as  shown  schematically  in  Figure  17.  The  shroud  Is 
being  designed  to  provide  a  sonic  flow  at  the  location  of  the  base  of  the  model, 
with  a  high  subsonic  flow  at  the  cavity  location.  An  analysis  is  underway  to  deter¬ 
mine  the  effect  of  heat  loss  from  the  flow  to  the  model  and  shroud  surfaces.  The 
mold  design  for  the  foams  for  these  models  is  in  preparation.  The  tests  will  be 
run  for  from  5  to  10  seconds  to  obtain  approximately  one  inch  of  recession  of  the 
foam  surface. 

2.  Phase  tl_ 

A  review  of  available  high-temperature  lightweight  foams  is  in  progress. 
The  seminar  on  Cellular  Plastics  was  attended  at  the  University  of  Michigan 
(August  3-7).  As  a  result  of  this  and  other  contacts,  additional  new  foams  have 
been  added  to  the  list  of  promising  candidates.  Samples  of  a  thermally  stable 
2  lb/ft^  high-temperature  urethane  foam  have  been  received  from  the  Amoco 
Chemical  Corporation.  A  list  of  foam  candidates  for  the  screening  tests  is  given 
in  Table  10. 

The  fabrication  procedures  for  producing  urethane  foams  of  CPR-1B  have 
been  determined,  and  density  variations  have  been  established.  By  removing  the 
outer  one  inch  areas  of  a  foamed  block,  the  variations  in  density  are  held  to  within 
0.  5  lb/ft^.  This  foam  will  be  supplied  to  the  Naval  Research  Laboratory  in 
12-inch  cubes  for  hypervelocity  impact  testing. 
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The  Linde  torch  equipment  is  being  readied  for  the  screening  phase  of  the 
program.  Agrei  ment  was  reached  on  the  design  of  the  sample  holder  and  tempera¬ 
ture  monitoring.  Design  of  a  fixture  for  installing  thermocouples  in  the  foam 
samples  prior  to  exposure  to  the  Linde  torch  is  being  pursued  so  thRt  the  resultant 
data  may  be  meaningful  and  consistent.  This  fixture  should  be  able  to  double  as 
a  device  capable  of  measuring  recession  rates  in  future  testing,  Tests  will  be 
performed  to  characterize  promising  foams  by  determining  significant  temperature 
ranges  for  mass  loss  and  phase  transitions. 


TA13LF.  10.  FOAM  CANDIDATES  FOR  SCREENING  TESTS 


Class 

Specific  Material 

Density  lb/ft3 

Max  Temp.  IJF 

1. 

Epoxy 

Ecco-foam  EFB 

10,  15 

300 

2. 

Epoxy 

Ecco-foam  EFF 

10,  15 

300 

3. 

Silicone 

DCR-7002 

14 

--•500 

4. 

Ceramic 

Ecco-foam  LM-43A 

15 

1000 

5. 

Ceramic 

Ecco-foam  Q  (or  equivalent) 

12 

2500 

6. 

Ceramic 

Foamglass 

9 

1000 

7. 

Urethane 

CPR  18 

8-15 

350 

8. 

Urethane 

CPR  18  +  Fire  Retard. 

8-15 

350 

9. 

Urethane 

CPR  18  +  Fire  Retard. 

+  Hollow  Glass  Fibers 

10-15 

350 

10. 

Urethane 

CPR  18  ♦  Fire  Retard. 

+  Flake  Glass 

10-15 

350 

11. 

Urethane 

CPR  18  +  Fire  Retard, 
i  Ecco-spheres  Si 

10-15 

350 

12. 

Urethane 

Hetrofoam  368/209 

4-15 

350 

13. 

Urethane 

Dalsawood  (Fire  Retardant 
Impregnant,  as  is  and  ns 
plywood) 

6-8 

14. 

Composite 

4  or  5  backed  up  with  15 

10-15 

>350 

15. 

Urethane 

Hetrofoam  HLR-250  •  PAPI 

4-15 

>350 

10. 

Composite 

Glass  Phenolic  Honeycomb 

4  Glass  backing,  all  backed 
with  Urethane  foam  15  or  7. 

■350 

17. 

Urethane 

Hetrofoam  HLR-250  (or 
equivalent)  +  PAPI  *  Hollow 
Glass  Fibers 

10-15 

>350 

18. 

Phenolic 

Tetrafonm 

10 

'350 

19. 

Urethane 

Amoco  Chemical  Corp. 

2 

500 
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Figure  18.  Arc  Shroud  Calorimeter  Model 
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m.  THERMOSTRUCTURAL  KILL  STUDY 

A  one-yeur  study  was  begun  in  August  with  the  objectives  of: 

1.  the  formulation  of  a  simplified  computational  model  for  the  determination 
of  thermostructural  vehicle  kill  due  to  hypervelocity  impact  perforations;  and 

2.  the  application  of  the  computational  model  to  selected  advanced  ICBM 
re-entry  vehicle  designs  of  two  general  classes:  (a)  unhardened;  and 

(b)  hardened. 


A.  Program  Plan 

This  aspect  of  the  continuing  work  under  the  HKM  program  comprises 
the  following  subtasks. 

Subtask  1:  Approximate  heat  shield  and  structure  designs  will  be  prepared  for 
several  advanced  ICBM  re-entry  vehicle  configurations  having  ballistic  co¬ 
efficients  greater  than  1500  psf.  Cone  angle  and  length  will  be  varied  over 
representative  ranges  of  interest.  Three  types  of  structure  design  (e.  g. ,  monocoque, 
ring- stiffened,  and  honeycomb-sandwich)  and  two  materials  (e.  g. ,  steel  and 
aluminum)  will  be  considered  for  each  configuration.  The  designs  will  be  "un¬ 
hardened"  per  the  following  definition:  no  foam  within  the  internal  compartments 
of  the  vehicle,  standard  aft  cover  design,  and  no  nuclear  blast  loading  requirement. 

Subtask  2:  Approximate  structure  designs  will  be  prepared  for  selected 
configurations  of  Subtask  1  which  are  "hardened"  per  the  following  definition: 
foam  within  the  Internal  compartments  of  the  vehicle  for  shell  stiffening,  aft 
cover  designed  to  withstand  cone  pressure,  and  various  levels  of  nuclear  blast 
loading  requirement. 

Subtask  3:  Structural  failure  criteria  will  be  formulated  for  forecone  and  aftcone 
regions  which  have  been  perforated,  for  unhardened  and  hardened  designs.  These 
criteria  will  be  based  primarily  on  existing  analytical  techniques,  with  a  minimum 
of  additional  cylinder  buckling  tests  similar  to  those  reported  on  in  Reference  1 
being  performed  if  required. 

SubtaBk  4:  An  approximate  model  will  be  formulated  for  computing  structure 
temperature  rises  and  foam  ablation  rates  due  to  internal  heating  of  perforated 
compartments.  The  model  will  include  the  effects  of  perforation  growth  during 
re-entry  and  non-uniform  internal  heating  distribution.  Simplified  relations  will 
be  employed  in  order  to  express  the  internal  heating  rates  as  functions  of  free- 
stream  quantities.  Approximate  trajectories  and  pressure  distributions  will  be 
used  in  an  attempt  to  obtain  closed  form  solutions  for  Internal  pressure  and 
structural  temperature  rise  as  functions  of  ballistic  coefficient,  re-entry  velocity 
and  path  angle,  and  altitude  at  which  perforation  occurs. 

Subtask  5:  Parametric  computations  will  be  made  for  the  vehicle  designs  of 
Subtasks  1  and  2  to  determine  perforation  size  required  for  thermostructural  kill 
as  a  function  of  altitude  at  which  perforation  occurs  and  altitude  of  kill.  The 
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computations  will  assume  a  single  perforation  in  either  the  forecone  or  aftcone 
region.  Approximate  heating  loads  to  exposed  payload  and  components  following 
thermostructural  kill  will  also  be  determined.  Results  for  the  "unhardened"  class 
of  vehicle  will  include  the  effect  of  5  and  10  percent  weight  additions  in  the  form 
of  foam  for  thermal  hardening  only. 


B.  Progress  to  Date 


Subtask  1 


Initial  study  under  this  task  has  resulted  in  the  selection  of  seven  vehicle 
configurations  having  ballistic  coefficients  from  1700  to  4400  psf,  with  cone  half¬ 
angles  from  8  to  12  degrees  and  lengths  from  100  to  267  inches.  A  summary  of 
these  configurations  is  given  in  Table  11. 

TABLE  11.  SUMMARY  OF  VEHICLE  CONFIGURATIONS 


Configuration 

Number 

Cone 

Half-Angle 

(degrees) 

Length 

(Inches) 

w/cda 

(psf) 

1 

8 

125.0 

2470 

2 

8 

153.5 

4375 

3 

8 

267.3 

3410 

4 

10 

99.7 

1720 

5 

10 

133.7 

2675 

6 

10 

213.2 

2335 

7 

12 

176.9 

1690 

Subtaak  4 


A  review  of  calorimeter  heat  transfer  and  perforation  ablation  measure¬ 
ments  obtained  during  the  AEDC  Tunnel  C  and  the  Malta  Rocket  Exhaust  Internal 
Heating  Test  Programs  is  in  progress.  Turbulent  shear  layer  theory  will  be 
applied  to  this  data  to  formulate  a  correlating  parameter  for  the  rate  of  perforation 
enlargement  due  to  re-entry  heating. 
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ABSTRACT 


This  report  discusses  work  performed  on  contract  NONR  3307(00)  in 
thu  period  March  to  September  1964.  The  major  portion  of  the  effort 
expended  to  date  has  been  in  the  area  of  instrumentation  and  investi¬ 
gation  of  arc  instabilities.  Descriptions  are  presented  of  the 
apparatus  constructed  for  measurement  of  gas  enthalpy,  composition, 
velocity  and  density  of  a  free  compressible,  turbulent,  high  enthalpy 
jet. 
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NOMENCLATURE 


The  list  below  defines  the  principal  notations  used  in  the  report. 
Those  symbols  that  were  used  infrequently,  are  defined  in  the  text. 


Symbols 

CP 

E 

h 

1 

K 

K<ji 

M 

P 

r 

R 

T 

U 

W 

X 

is 

p 

Subscripts 

A 

C 

CI. 

g 


Specific  heat  at  constant  pressure 

Voltage  applied  to  arc 

Enthalpy 

Current 

Turbulent  Exchange  Parameter  for  Velocity 

Turbulent  Exchange  Parameter  for  Temperature 

Mach  number 

Pressure 

Radial  distance 

Jet  nozzle  radius 

Temperature 

Velocity 

Weight  flow  rate 
Axial  distance 
Specific  heat  ratio 
Density 

Subscripts 


Ambient  conditions 
Coolant 

Jet  centerline 
Gas  sample 


S  Stagnation  conditions 

0  Jet  exit  conditions 

1  Probe  entrance 

2  Probe  exit 
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The  work  to  be  accomplished  at  Avco/RAD  under  this  program  consists 
of  an  extensive  experimental  study  of  the  characteristics  of  a  high 
enthalpy  jet  mixing  with  a  stationary  ambient  atmosphere.  Velocity, 
enthalpy,  and  density  distributions  obtained  under  a  variety  of  jet 
operating  conditions,  along  with  pressure  and  heat  flux  distributions 
measured  on  a  flat  plate  held  normal  to  the  jet,  will  be  correlated 
with  existing  theories  in  order  that  a  capability  for  the  prediction 
of  heat  transfer  may  result. 

A  high  enthalpy  jet  is  produced  by  an  electric  arc  heater.  This 
jet  is  probed  with  an  enthalpy  probe,  which,  together  with  an  assumption 
that  the  gas  in  the  jet  is  in  thermal  equilibrium,  results  m  sufficient 
information  to  deduce  the  desired  enthalpy,  velocity,  and  density  values. 
The  resulting  data  on  turbulent  jet  mixing  is  intended  to  produce  values 
of  the  turbulent  exchange  coefficient  (related  to  the  eddy  viscosity)  in 
a  regime  of  temperature  and  Mach  Number  for  which  such  information  is 
not  presently  available. 

In  addition  to  the  determination  of  the  aerodynamic  properties  of 
the  flow  and  our  ability  to  predict  them,  the  heat  transfer  to  a  flat 
plate  placed  normal  to  the  flow  will  be  investigated.  Pressure  distri¬ 
bution  on  a  flat  plate  in  the  jet  will  be  measured  and  used  to  determine 
the  stagnation  point  velocity  gradients,  which  will  in  turn  be  u9ed  to 
calculate  the  expected  heat  flux.  These  calculated  values  of  heat  flux 
will  be  compared  with  measured  values. 

Ultimately  this  effort  is  intended  to  result  in  a  procedure  whereby 
the  heat  transfer  from  a  jet  of  specified  properties  to  a  flat  plate  held 
normal  to  the  jet  axis  may  be  predicted.  The  program  being  undertaken 
is  thus  primarily  experimental,  including  an  attempt  to  correlate  data 
with  prediction. 
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The  facility  employed  in  producing  iho  high  enthalpy, 
compressible  jet  to  be  investigated  i  the  1.5  Megawatt  Plasma  Generator 
shown  in  Figures  1  &  2.  Basically  the  irr  heater  consists  of  an  exit 
nozzle  anode  assembly,  on  ore  chamber,  and  a  cathode  assembly.  Although 
the  arc  can  be  modified  to  give  three  different  configurations  ill ,  the 
system  to  be  used  in  this  study  is  the  gas  stabilized  arc  configuration 
(with  the  option  of  magnetically  driven  arc  roots;  shown  schematically 
in  Figure  3.  The  facility  has  the  capability  of  being  operated  over  a 
large  range  of  chamber  pressures  (1)  thus  , malting  it  possible  to  produce 
both  subsonic  and  sonic  jets  with  various  degrees  of  undcrcxpansion .  At 
pressure  levels  to  be  considered  in  this  program  it  is  not.  necessary  to 
drive  the  arc  roots  with  magnetic  fields  hence  for  its  present  appli¬ 
cation  the  facility  may  be  considered  to  be  gas  stabilized. 

Gas  enters  the  plasma  generator  through  a  series  of  inlet 
holes  in  the  arc  chamber  and  is  heated  to  a  high  enthalpy  by  flowing 
through  the  discharge  which  passes  across  the  gap  between  the  cathode 
and  anode.  The  heated  gas  then  exhausts  to  the  atmosphere!  through  the 
exit  nozzle. 


The  working  fluid  to  be  employed  in  all  experiments  con¬ 
ducted  in  the  program  is  nitrogen.  The  use  of  air  was  not  considered 
since  it  would  result  in  erosion  of  the  cathode  assembly  thus  contamiat- 
,i  the  exhaust  jet.  Simulated  air  (i.e.,  using  nitrogen  as  the  primary 
fluid,  passing  it  through  the  arc  column,  and  then  adding  oxygen  to  the 
hot  gas  in. a  plenum  chamber  between  the  anode  and  exit  nozzle  to  produce 
simulated  air)  was  not  considered  since  the  amount  of  mixing  using  t.his 
technique  has  not  been  established;  hence  it  would  not  be  possible  to 
insure  that  the  exhaust  jet  would  be  homogeneous  in  nature-  Nitrogen  as 
a  working  fluid  has  several  advantages  in  that  the  arc  may  be  run  for 
extended  periods  of  time,  the  amount  of  contamination  is  kept  to  a  mini¬ 
mum,  and  it  simulates  air  to  a  reasonable  degree. 

A  major  portion  of  the  investigations  in  this  program,  i.e., 
measurement  of  axial  and  radial  profiles  of  jet  density  velocity,  com¬ 
position  and  enthalpy,  are  to  bo  conducted  with  an  0  75-inch  diameter 
exit  nozzle.  Hence,  the  jet  will  be  considerably  larger  than  the  diag¬ 
nostic  probes  (described  below)  required  to  obtain  these  quantities. 
Operation  of  the  arc  facility  takes  place  in  a  region  of  mass  flow  rates 
and  gas  enthalpies  such  that  the  jet  Reynolds  Number  is  well  above  the 
value  of  approximately  1500  (2)  required  for  production  cf  a  turbulent 
jet . 

Operating  conditions  for  the  facility  are  tabulated 

be  low. 
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Figure  2 — Closeup  View  of  1.5  Megawatt  Plasira  Generator  Electrodes 


TABLE  I 


1.5  MEGAWATT  PLASMA  GENERATOR  OPERATING  CHARACTERISTICS 


Voltage,  volts 
Current,  amperes 
Powor,  kilowatts 
Efficiency,  % 


500-750 

500-1500 

250-1125 

45-60 


EXIT  NOZZLE  DIAMETER,  0.75-INCH 


Gas  Enthalpy,  H/RT0 
Gas  Temperature,  °K 
Gas  Density,  lb/ft3 
Mach  No. 

Gas  Velocity,  ft/sec. 
Reynolds  Number 
Exit  Pressure  Ratio 


75-165 
4000-7000 
5  X  10-3-2.5  X  10 
.9  -  1.0 
3750-5000 
18,000-8,000 
1-2 
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In  the  determination  of  mixing  parameters  which  would 
allow  one  to  predict  jet  decay  it  is  extremely  important  that  mix¬ 
ing  of  the  jet  and  ambient  fluid  should  be  the  result  of  turbulent 
shear  mechanisms  and  not  caused  by  jet  instabilities  peculiar  to  the 
particular  device  employed.  Relatively  large  scale  jet  fluctuations 
such  as  flapping  and/or  "on-off"  instabilities  may  cause  a  probe 
mounted  in  the  stream  to  be  periodically  immersed  m  the  jet  and 
ambient  fluids.  Due  to  the  response  time  of  the  probes,  the  measured 
quantities,  such  as  impact  pressure  gas  enthalpy,  and  composition, 
would  be  time  averages  of  the  jet  and  arr.bient  fluid  variables.  Data 
of  this  nature  would  indicate  that  the  axial  decay,  and  radial  spread 
of  the  jet  would  be  more  rapid  than  that  of  its  steady  counterpart. 

Hence,  mixing  parameters  obtained  with  an  unsteady  jet  would  be  ex¬ 
tremely  difficult  to  interpret.  Therefore-  it  is  essential  that  the 
jet  have  a  high  Strouhal  number,  i.e.,  any  instabilities  in  the  flow 
should  have  a  low  amplitude  and  a  high  frequency. 

An  investigation  of  the  time  dependent  characteristics  cf 
the  1.5  Megawatt  Plasma  Generator  exhaust  jet  has  been  initiated. 
Stability  of  the  high  enthalpy  jet  was  examined  by  high  speed,  motion 
picture  photographs  (2,000,  4,000,  and  8,000  frames/second).  The 
initial  studies  were  performed  with  the  arc  having  the  configuration 
shown  in  Figure  3.  A  typical  sequence  of  film  (8,000  frames/second) 
from  this  test  series  is  presented  in  Figure  4.  The  films  demonstrated 
that  the  jet  had  both  flapping  and  on-off  instabilities.,  both  of  which 
were  believed  to  be  associated  with  the  location  of  the  anode  arc  roots. 

A  significant  portion  of  the  arc  discharge  was  occurring  within  the 
nozzle  throat  and  in  the  vicinity  of  the  downstream  face  of  the  exit 
nozzle,  This  phenomena  was  responsible  for  the  bright  region  immediately 
downstream  of  the  nozzle  (see  Figure  4  i . 

Subsequent  to  these  experiments,  the  arc  configuration  (as 
shown  in  Figure  3.1  was  modified  by  the  addition  of  a  plenum  chamber 
between  the  anode  and  nozzle.  In  addition,  the  length  of  the  exit 
nozzle  was  increased.  Both  alterations  made  it  possible  to  eliminate 
the  discharge  occuring  within  the  throat  and  to  move  the  anode  arc  roots 
upstream  of  the  plenum  chamber.  The  location  of  the  plenum  between  the 
arc  discharge  and  exit  nozzle  made  it  possible  to  reduce  the  jet  in¬ 
stability.  A  typical  sequence  of  high  speed  motion  picture  film  (8,000 
frames/second)  for  the  modified  generator  is  included  in  Figure  4.  Al¬ 
though  there  are  some  fluctuations  still  present  in  the  jett!' 
flapping  instability  has  been  essentially  eliminated  and  "on-off"  in¬ 
stabilities  have  been  considerably  reduced. 

At  present  the  arc  is  being  further  modified  by  addition 
of  a  weak  magnetic  field.  It  is  expected  that  the  field  will  help  to 
produce  a  more  diffuse  arc  column  and  honce  increase  the' j et  stability. 

Further  information  on  the  mechanism  of  jet  decay  and 
stability  can  be  obtained  by  measurement  of  the  jet  acoustic  power 
spectrum.  Jet  mixing  caused  by  turbulent  shear  mechanism  produces  a 
"white"  noise,  i.e.,  the  sound  pressure  level  (SPL)  would  be  constant 
over  a  wide  range  of  frequencies,  whereas  mixing  produced  by  jot 
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Figure  *i«-  High  Speer]  Film  Sequence  of  1.':'.  MW  Plasma  Generator  Exhaust  Jet. 

Without  Plenum  (left),  With  Plenum  (right):  Camera  Gpee.l  (i-'.OOO  fpr) 
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instability  results  in  a  peak  in  the  SPL  at  various  frequencies  associated 
with  the  particular  apparatus.  ; 

Measurements  were  made  to  determine  the  SPL  of  the  jet 
produced  bv  the  1.5  MW  Arc  facility  (with  plenum  chamber).  In  addition 
to  the  alterations  described  above  the  arc  was  modified  further  by 
addition  of  a  large  diameter  plate  in  the  exit  plane  of  the  nozzle  so 
as  to  simulate  a  jet  issuing  from  a  hole  in  an  infinite  wall.  Measure¬ 
ments  were  made  at  arc  operating  conditions  identical  to  those  employed 
in  the  last  series  of  high  speed  motion  picture  tests  (see  Figure  4.). 

Data  obtained  were  limited  to  tape  recordings  of  the  jet  noise  which 
could  then  ho  analyzed  in  the  laboratory  The  following  equipment  was 
employed  in  these  experiments t 

Ampex  Tape  Recorder,  Model  601 

Massa  Microphone,  Model  M-213,  S/N  33 

Massa  Amplifier  and  Power  Supply,  Model  M-185 

Massa  Pre-Amplifier.  Model  M-114B,  S/N  572 

Massa  Microphone  Base,  Model  M-218 

The  microphone  was  oriented  for  normal  incidence  with  respect  to  the  jet 
nozzle  and  in  the  plane  of  the  baffle  plate.  The  microphone  was  clamped 
on  a  stand  which  rested  on  the  concrete  floor.  Sheets  of  rubberized 
horse-hair  packing  material,  approximately  two  inches  thich,  were  placed 
on  the  floor  in  the  vicinity  of  the  microphone  so  as  to  reduce  sound 
reflection  from  the  floor  to  the  sensor.  Facing  upstream,  the  microphone 
was  located  in  the  third  quadrant. 

Recordings  obtained  in  each  test  were  analyzed  with  a  1/3 
octave  frequency  analyzer  by  running  a  continuous  loop  of  the  tape  through 
the  spectrum  analyzer  and  examining  the  noise  level  obtained  with  each 
filter  separately.  Two  loops  were  analyzed  for  each  experiment  and  the 
results  were  corrected  to  that  which  would  be  obtained  with  constant  band¬ 
width  filters.  The  output  of  a  General  Radio  Noise  Generator  Model  1390B 
was  recorded  and  analyzed  to  determine  the  deviation  of  the  tape  recorder 
from  flat  response.  Noise  level  of  all  equipment  associated  with  the  arc 
facility  (without  the  arc  running)  was  also  measured  and  analyzed. 

The  results  of  these  experiments  in  the  form  of  SPL  or  db 
with  reference  to  0.002  dyne/cm2  as  a  function  of  frequency  are  presented 
in  figures  5  through  7.  Comparison  of  the  data  obtained  on  the  noise 
level  of  the  auxilliary  arc  equipment  and  that  produced  by  the  jet  indicate 
that  the  rise  in  the  SPL  of  the  jet  at  frequencies  below  100  cps  and  above 
8000  cps  may  be  the  result  of  ambient  noise.  The  results  indicate  that 
the  jet  audio  output  was  relatively  flat  over  a  wide  range  of  frequency. 

No  single  bandwidth  had  noise  levels  greatly  different  from  the  mean  SPL 
indicating  that  no  resonant  characteristic  of  the  arc  showed  up  significantly 
in  the  exhaust  jet. 
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B  GAB  PROBES 


For  a  considerable  period  of  time  a  company  sponsored  probe  develop- 
rent  proqram  has  been  underway  at  Avco/RAD.  The  purpose  of  this  program 
has  been  to  develop  thermodynamic  and  calorimetric  probes  for  diagnostic 
studies  in  all  facilities  which  provide  high  enthalpy  gases  for  laboratory 
use.  Probes  have  been  developed  to  determine  local  enthalpies,  gas  com¬ 
positions,  mass  flow  rates,  heat  fluxes,  and  both  static  and  impact 
pressures . 

Of  the  facilities  available  for  test  purposes,  the  most  severe  conaition 
are  encountered  in  the  10  Megawatt  Arc(3).  For  that  reason  it  was  felt 
that,  if  the  probes  could  withstand  the  high  heat  flux  at  the  1U  Megawatt 
Arc  without  melting,  then  there  would  be  nc  questions  about  their  ability 
to  withstand  the  environments  produced  by  other  plasma  generators. 

Early  in  *-he  probe  development  program  it  was  decided  to  use  copper 
for  all  sir  races  that  are  in  a  contact  with  the  hot  gases.  Despite  the 
fact  that  copper  has  a  relatively  low  melting  point,  it  was  felt  that 
with  a  proper  design  of  a  water-cooling  system,  advantage  could  be  taken 
of  the  high  thermal  conductivity  of  this  metal. 

1 . i  Prototype  Water  cooled  Probe 

In  the  initial  phase  of  this  probe  development  program  a  proto¬ 
type  probe  (Figures  8  &  B)was  constructed  and  placed  in  the  10  Megawatt 
Arc  exhaust  jet. Experiments  performed  with  an  0.25-inch  diameter  probe 
proved  successful  in  that  the  device  could  withstand  estimated  heat 
fiuxes  of  18  X  106  Btu/ft2hr  with  a  water  coolant  pressure  of  280  psig. 

In  view  of  the  favorable  experience  in  the  10  Megawatt  facility 
a  series  of  watercooled  calorimeters  for  measuring  heat  flux  have  been 
designed  and  are  in  operation.  The  calorimeters  vary  in  size  from  h  in. 
diameter  and  Larger-  depending  on  the  requirement.  Emphasis  is  being 
f..aced  in  future  design  cn  miniaturization,  fast  response  and  reduction 
f  error  due  to  heat  losses 

2.)  Mass  Flow  Probes 


Probes  designed  to  measure  the  local  mass  flux  in  a  supersonic 
gaa  ir-  have  been,  gpnstructed  and  tested  in  the  exhaust  of  a  30 

Kilowatt  arcjet  engine.  L*'Ihe  mass-flux  probes  were  0.250  inch  in  cuter 
diameter  and  had  various  diameter  (0.042-0.062  inch)  gas  sampling  tubes 
located  on  the  probe  axis.  In  contrast  with  the  prototype  calorimetric 
probe,  the  mass  flux  probes  have  a  conical  front  face  which  comes  to  i 
sharp  edge  at  the  junction  with  the  gas  sampling  tube.  Cone  half  angles 
used  in  these  probes  were  20°  and  45°.  A  typical  mass  flux  probe  is 
shown  in  Figure  10, 

Experiments  conducted  with  the  masB  flux  probe  (200  nose  angle1' 
were  performed  in  the  30  kilowatt  arcjet  engine.  The  probe  was  plated 
as  close  as  possible  to  the  exit  plane  of  the  engine  and  operated  success- 
fully  to  within  0.10  inch  of  the  jet's  axis  of  symmetry.  Experience  with 
this  probe  indicated  that  when  a  reasonable  pressure  difference  was  main¬ 
tained  across  the  probe  sampling  tube,  all  shocks  in  the  stream  tube 
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Figure  H — Prototype  Water  Cooled  Probe 


Probe,  Supersonic  FLov 


being  sampled  occur  within  the  probe  i.e.,  the  probe  shock  was 
"swa 1 lowed  ' . 

A  typical  mass  flux  profile  obtained  in  the  arcjet  engine  exhaust 
by  using  water  cocied  and  carbon  mass  flux  probes  in  a  plane  0,0625  inch 
from  the  nozzle  exit  is  presented  in  Figure  11.  The  nature  of  the  mass 
flux  profile  was  in  agreement  with  prediction  and  the  total  integrated 
mass  flux  was  within  17%  of  the  flow  to  the  engine. 

Future  experiments  to  be  conducted  with  the  supersonic  mass 
flux  probes  are  to  investigate  the  effects  of  probe  nose  angle  and 
sampling  tube  diameter  on  the  ability  to  swallow  the  probe  bow  shock. 

3- )  Pressure  Probes 


Flat  nosed  impact  pressure  probes  were  fabricated  and  tested  in 
the  10  Megawatt  facility  exit  jet,  In  general  the  pressure  probe  was 
of  the  same  shape  as  the  prototype  water-cooled  probe  (Figures  8  &  9) . 

The  pressure  probe  operated  satisfactorily  at  heat  flux  levels  up  to 
18  X  106  BtU./ft2hr. 

In  addition  a  water  cooled  probe  (Figure  12)  designed  to  measure 
static  pressure  within  a  high  enthalpy,  high  velocity,  gas  stream  has 
been  developed  In  brief,  the  probe  has  four  pressure  ports  on  the  side 
of  the  instrument  to  measure  3tatic  pressure  within  a  plane  perpendicular 
to  the  direction  of  gas  flow.  The  openings  are  equally  spaced  around  the 
surface  of  the  probe.  Pressure  lines  leading  from  each  pressure  port 
open  into  a  common  plenum  within  the  probe  body,  hence  the  static  pressure 
measured  with  this  device  wil .  be  an  average  value. 

4  ;  Enthalpy  Probes 

The  neasurenent  of  gas  properties  at  atmospheric  pressures  and 
very  high  temperatures  has  long  been  a  difficult  problem  Recently 
water  cooled  probes  capable  of  steady  state  operation  in  such  environ¬ 
ments  have  been  developed  '.i)  These  probes  have  the  capability  of 
determining  total  pressure,  composition  and.  m  some  cases,  mass  flux 
as  well  ae  enthalpy  A  diagram  of  Buch  a  probe  is  given  in  Figure  13 
Such  probes  are  based  on  the  principle  that  the  rate  of  heat  loss  by  the 
gas  flowing  through  the  probe  is  eque ’  to  the  rate  of  heat  input  to  the 
cooling  water  Errors  due  to  heat  l  t  from  gas  which  does  not  flow 
through  the  probe  are  effectively  eliminated  by  performing  a  '‘tare’' 
measurement  that  is,  by  recording  the  coolant  temperature  rise  when  no 
gas  flows  through  the  probe  and  subtracting  it  from  the  temperature 
rise  cf  the  coolant,  when  the  sample  is  allowed  to  pass. 

Two  types  of  enthalpy  probe  have  received  primary  consideration 
for  use  in  this  progran. .  The  first  is  termed  the  "single- jacketed " 
enthalpy  probe.  A  probe  of  this  type  was  developed  by  J  Grey  at  Princeton 
and  is  described  in  detail  in  Reference  5  One  of  these  probes  was  pur¬ 
chased  from  the  Greyrad  Corp  In  addition  a  probe  of  this  type  was  con¬ 
structed  at  A vco  in  order  to  allow  the  use  of  lower  water  pressure  than 
are  required  for  operation  cf  the  Greyrad  probe. 
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Figure  11 — Mass  Flux  Profile,  Avco  Plasma  Ergine 
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Fitn.ire  12 — Gt.at.ie  Pressure  Probe 


The  Avco  single  jacketed  enthalpy,  Figure  14,  0.187  inch  o.d., 
was  tested  in  the  exhaust  of  the  30  kilowatt  arcjet  engine.  The  probe 
was  initially  tested  six  inches  downstream  from  the  engine  orifice  and 
then  moved  gradually  to  three  inches.  The  impact  pressure  at  the  probe 
tip  was  sensed  without  difficulty,  and  the  coolant  temperature  rise  was 
also  measured  with  sufficient  accuracy  once  the  coolant  flow  rate  was 
reduced  to  a  low  enough  value.  This  initial  experiment  with  the  single- 
jacketed  enthalpy  probe  proved  that  reasonable  values  of  stagnation 
enthalpy  may  be  determined  with  this  device. 

Nevertheless,  in  all  runs  made  with  this  probe,  the  tare  reading, 
or  coolant  temperature  rise  with  no  flow  through  the  probe,  was  a  mini¬ 
mum  of  80%  of  the  coolant  temperature  rise  with  the  probe  opened.  Be¬ 
cause  of  the  necessity  of  subtracting  large  numbers  during  the  enthalpy 
calculation,  it  is  difficult  to  visualize  the  possibility  of  accurate 
determinations  of  the  total  enthalpy  with  this  device.  It  is  obvious 
that  a  small  error  in  measuring  the  tare  will  cause  a  relatively  large 
error  in  computing  the  enthalpy.  Such  small  errors  in  determining  the 
tare  may  be  inherent  in  the  technique  since  the  tare  which  is  measured 
when  there  is  no  gas  flow  through  the  tube  is  not  the  same  "tare"  that 
exists  when  there  is  flow.  The  difference  between  these  two  tares, 
however,  may  be  so  small  that  for  practical  purposes  it  may  be  ignored. 

In  order  to  reduce  these  difficulties  as  much  as  possible  and 
to  gain  the  greatest  accuracy  possible  in  the  enthalpy  determination,  a 
double-jacketed  probe  termed  a  "split-flow"  enthalpy  probe  was  designed 
and  constructed.  This  probe  is  shown  in  Figure  15.  This  probe  is  0.218 
inch  o.d.  and  is  believed  to  incorporate  a  significant  improvement  over 
the  single-jacketed  probe.  In  this  probe  the  tare  values  can  be  reduced 
to  an  insignificant  portion  of  the  total  heat  gain  by  the  water  and  thus 
greater  accuracies  may  be  expected  in  determining  the  enthalpy  of  the 
gas.  A  schematic  drawing  of  this  Avco  probe  is  shown  in  Figure  16. 

The  probe  consists  of  four  concentric  tubes.  The  innermost  tube 
is  open  to  permit  gas  flow  through  it.  The  annulus  between  the  second 
and  third  tubes  is  used  to  admit  the  water  coolant.  At  the  tip  of  the 
probe  the  water  flow  is  split,  hence  the  designation  "split-flow".  The 
outer  water  is  used  to  cool  the  outer  probe  surface  while  the  inner 
water  is  utilized  to  cool  the  gas  which  is  tested.  With  this  design 
the  tare  is  proportional  only  to  the  temperature  rise  of  the  test-gas 
water  coolant  when  there  is  no  gas  flow  through  the  center  of  the  probe. 
Since  the  outer  water  flow  rate  can  be  adjusted  so  that  its  temperature 
rise  is  equal  to  that  of  the  inner  water,  there  should  be  no  heat  trans¬ 
ferred  between  the  outer  water  colant  and  the  test-gas  coolant. 

For  mass-flow  measurements,  the  tip  of  the  probe  can  be  made 
knife-edged  to  "swallow"  the  bow  shock  when  the  flow  is  supersonic. 

It  is  expected  that  the  same  probe  can  be  used  to  determine  impact 
pressures,  mass  flow  rates  and  enthalpies. 

Before  «n  enthalpy  probe  is  used,  it  should  be  calibrated.  When 
an  arc  facility  is  employed,  the  arc  itself  may  be  utlized  as  the  stand¬ 
ard  calorimeter.  In  a  instrumented  arc  generator  the  total  energy  to 
the  plasma  may  be  determined  by  subtracting  the  losses  from  the  total 
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electrical  input.  This  value  is  then  compared  to  the  integrated  value 
of  the  individual  points  indicated  by  the  probe  when  the  plasma  jet  is 
traversed  radially.  The  difference  between  the  integrated  value  and 
the  computed  input  to  the  plasma  is  assumed  to  be  the  probe  error. 

The  mathematical  relationshi  for  determining  the  percent  error  is 

L  (  i  aiT  r  pu  W  dir  )  -  C  EX  -  wcCp iqq 
EX  -  Cp  4>Tc 


as  follows: 
f-  r  re r 


The  relative  merit  of  the  ‘split -flow"  probe  as  compared  to 
the  "tare"  instrument  was  determined  by  testing  the  "split-flow"  device 
in  the  Model  500  ARC  facility  (6) .  This  facility  produces  a  subsonic 
exhaust  jet  at  gas  enthalpies  ranging  from  3,000  to  11,000  Btu/lb.  with 
stagnation  point  heat  transfer  rates  of  2,000-4,000  Btu/ft2  second  to 
conical  bodies  having  a  spherical  nose  radius  of  0.125  inch.  The  energy 
balance  technique  used  to  determine  the  exhaust  gas  enthalpy  has  been 
found  to  be  accurate  for  this  facility, (6)  hence  the  Model  500  is  believed 
to  be  an  excellent  device  for  probe  calibration. 

A  "split-flow"  probe  (0.218  in  outside  diameter)  was  placed  in 
the  Model  500  exhaust  Jet.  In  this  particular  series  of  tests,  nitrogen 
was  employed  as  the  working  fluid.  The  arc  heated  gases  exhausted  to 
the  atmosphere  through  a  nozzle  having  an  exit  diameter  of  inch  as 
a  subsonic  Jet  with  a  Mach  number  of  0.70.  A  series  of  enthalpy  measure¬ 
ments  were  made  at  axial  distances  between  9  and  24  nozzle  radii  (Figures 
17  and  18) .  The  probe  was  not  moved  closer  to  the  exit  nozzle  since 
the  measured  temperature  rise  of  the  probe  coolant  indicated  that  boil¬ 
ing  was  imminent  at  the  smaller  axial  distance.  This  heat  transfer 
limitation  of  the  probe  was  a  result  of  the  low  water  pressure  available 
at  the  facility  (200  psig)  and  the  restriction  of  the  coolant  flow  by  a 
0.006  inch  outer  water  passage  within  the  probe.  The  "split-flow"  probe 
has  been  redesigned  to  provide  greater  coolant  flow  through  the  outer 
water  channel.  Two  probes  of  this  type  have  been  constructed  and  will 
be  used  in  conjunction  with  a  high  pressure  water  supply  system  (200- 
1,000  psia)  for  the  remainder  of  this  program.  These  probes  are  to  be 
calibrated  using  the  Model  500  Arc  as  a  reference. 

It  is  significant  that  tare  values  considerably  less  than  the 
minimum  of  80%  found  with  the  single-jacketed  probe  were  recorded  during 
these  tests.  A  facsimile  of  a  typical  recorder  trace  obtained  during 
operation  of  the  split-flow  probe  in  the  arc  facility  is  shown  in 
Figure  19.  Tare  measurements  were  reduced  to  values  between  30  and  60% 
of  the  final  gas  coolant  temperature  rise.  This  was  so  even  though  no 
effort  was  made  at  that  time  to  minimize  the  tare  by  adjustment  of  the 
inner  and  outer  flow  rates.  Since  this  feature  is  expected  to  provide 
a  greater  increase  in  the  accuracy  of  the  enthalpy  determinations  at 
the  expense  of  very  little  added  probe  complexity,  the  Avco  "split-flow" 
enthalpy  probe  has  been  chosen  as  the  primary  experimental  tool  to  be 
employed  in  this  program. 

Although  the  positioning  of  the  probe  within  the  arc  exhaust 
jet  was  not  to  the  accuracy  required  for  the  measurements  to  be  made  in 
this  program,  jet  mixing  parameters  were  determined  using  the  method 
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by  Warren(7).  The  analytical  curves  (Figures  17  and  18)  were  found  by 
extrapolating  and  irterpolating  between  Warren's  results  for  the  appropriate 
temperature  ratio  and  Mach  number.  Values  of  the  mixing  parameters  K  and 
Kt  were  determined  to  provide  the  best  fit  to  the  Model  500  data.  Calculated 
mixing  parameters  found  were  a  factor  of  two  larger  than  predicted  with 
Warren's  correlation.  It  should  be  noted  that  since  the  accuracy  of  the 
positioning  device  used  in  these  tests  is  an  unknown  factor,  and  that  it 
was  not  possible  to  use  the  enthalpy  probe  close  to  the  jet  exit  plane 
the  results  obtained  in  the  Model  500  experiments  should  be  viewed  with 
caution. 

4. )  Calorimetric  Probes 


The  calorimetric  probes  are  used  to  determine  the  heat  flux  from 
a  hot-jet  which  comes  in  contact  with  the  calorimeter  surface.  From 
equationd  derived  from  theoretical  heat  transfer  analysis (8)  it  is  also 
possible  to  compute  the  enthalpy  of  the  gas  when  the  pressure  and  heat 
flux  are  known.  It  is  to  be  noted  that  whenever  these  calorimetric  probes 
are  used  in  a  reacting  gas,  the  catalytic  activity  of  the  wall  surface 
must  be  considered  (reference  9) .  Some  experimenters  believe  that  the 
the  discrepancies  in  certain  reported  data  are  due  to  the  differences  in 
the  catalytic  activity  of  the  surfaces  of  the  calorimetric  devices  used 
in  the  various  experiments. 

It  is  convenient  to  divide  these  calorimeters  into  two  main  groups: 
(a)  Steady  state  calorimeters  and  (b)  transient  state  calorimeters. 

(a)  Steady  State  Calorimeters 

The  steady  state  calorimeters  are  usually  made  of  copper 
and  are  water  cooled.  They  are  generally  designed  so  that  only  the  sur¬ 
face  normal  to  the  Jet  is  exposed  to  the  gas.  To  insure  this  condition, 
a  water  coo.led  copper  guard  is  used.  Ab  shown  in  Figure  20,  the  calori¬ 
meter  and  guard  are  in  the  form  of  two  concentric  tubes  closed  at  the 
ends . 


The  heat  flux  is  determined  from  measurements  of  the  water 
flow  rate,  temperature  difference  between  water  in  and  water  out  of  the 
calorimeter,  and  the  area  of  the  calorimeter  surface  exposed  to  the  he*-, 
gas.  (This  method  is  sometimes  called  the  "cold  wall  method"  because  f 
the  relatively  cool  contact  surface) .  The  temperature  of  the  water  in 
the  guard  is  measured  and  the  flow  rate  is  adjusted  to  minimize  the 
temperature  difference  between  the  adjacent  wall  of  the  calorimeter  and 
the  guard  ring. 

Calorimeters  of  this  type  have  been  used  at  this  laboratory 
for  several  years  with  relatively  good  results.  However,  it  is  believed 
that  they  must  be  re-dosigned  to  reduce  heat  loss  errors  and  also  miniatur¬ 
ized  so  that  local  values  may  be  determined.  The  re-design  of  this 
calorimeter  is  in  progress  at  the  present  time. 

Another  type  of  steady  state  calorimeter  is  now  being 
developed.  In  this  calorimeter,  as  shown  in  Figure  20,  a  resistor  to 
heat  transfer,  in  this  case  constantan , is  sandwiched  between  two  copper 
plates.  One  of  the  copper  plates  is  water  cooled  while  the  other  plate 
is  in  contact  with  the  hot  gas.  The  copper  plates  and  the  constantan  are 
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in  effect  a  coppcr-constantan  thermocouple  with  the  "hot"  junction  on 
the  gas  side  and  the  "cold"  junction  of  the  water  aide  of  the  eonstantan. 
By  connecting  the  two  copper  plates,  through  copper  leads,  to  a 
potentiometer  the  emf  developed  due  to  the  temperature  difference  at 
the  junctions  is  measured.  Since  the  emf  is  approximately  proportional 
to  the  neat  flux,  the  calorimeter  can  be  calibrated  to  indicate  the 
heat  flux  supplied  to  it. 

(b)  Transient  Calorimeters 


Transient  calorimeters  have  been  used  successfully,  at 
this  laboratory  for  quite  some  time  (Reference  10) .  The  heat  flux 
is  measured  py  determining  the  rate  ot  temperature  rise  of  the  surface 
of  the  solid  exposed  to  the  hot  gas.  In  general,  the  calorimeter 
consists  of  a  copper  cylinder  of  known  mass  and  dimensions.  One  end 
surface  of  the  cylinder  is  exposed  to  the  hot  gas  while  the  other  sur¬ 
faces  are  insulated.  It  is  aeeumed  that  the  heat  flow  is  unidirectional 
and  that  all  unexposed  surfaces  are  adiabatic.  The  disadvantage  of  this 
calorimeter  is  that  it  can  be  used  only  for  a  brief  period  of  time  (up 
to  several  seconds)  bei'.re  corrections  must  be  made,  in  some  cases,  for 
heat  losses  at  the  calorimeter  rear  face.  Its  important  advantage  is 
that  it  can  be  readily  used,  because  of  its  simple  design,  on  various 
shapes  to  measure  local  heat  transfer  rates. 

Based  on  the  above  described  calorimeters,  the  determination 
of  heat  fluxes  can  be  readily  made  and  cross-checked  using  the  steady 
state  and  transient  sidle  calorimetric  techniques. 

Probes  to  be  employed  specifically  on  this  program  include 

the  following: 

a. )  Split  Flow  Enthalpy  Probes  to  determine  the  axial  and  radial 
variations  in  jet  properties,  i.e.,  enthalpy,  velocity,  composition, 
and  density. 

b. )  Static  Pressure  Probes  to  measure  the  jet  static  pressure  which 
will  be  used  in  conjunction  with  the  impact  pressure  measurements  tc 
determine  local  Jet  velocity. 

c. )  Calorimetric  probes  (transient  device)  to  measure  the  heat 
flux  distribution  obtained  from  high  enthalpy  gas  jet  impinging 
upon  a  plate.  In  addition,  the  heat  flux  as  determined  by  steady 
water  cooled  probes  on  the  jet  axis  of  symmetry  can  be  employed 
with  standard  heat  flux  correlations (8)  to  determine  gas  enthalpy. 

The  latter  measurements  would  be  used  to  provide  a  check  on  enthalpy 
determinations  made  with  the  split  flow  enthalpy  ves. 

C.  Probe  Support  Equltment 

In  order  to  use  the  probes  described  previously  to  their  full  ad¬ 
vantage  various  pieces  of  support  equipnent  such  as  a  probe  position¬ 
ing  mechanism,  a  high  pressure  water  collant  system,  and  a  gas  collect¬ 
ing  system  have  been  constructed.  A  brief  description  of  these  items 
is  given  below. 
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A  .  )  Probe  Positioning  Mechanism 


In  order  to  obtain  reliable  and  reproducible  probe  measure¬ 
ments  it  is  important  to  have  a  means  of  positioning  the  probe  in  some 
known  location  relative  to  certain  convenient  references  such  as  the 
arc  nozzle  exit  plane  and  nozzle  axis.  To  accomplish  this  end  a  probe 
positioning  table  has  been  constructed, *a  photograph  of  the  table  is 
shown  in  figure  2.  With  this  table  he  various  probes  can  be  moved  in 
all  three  directions  in  a  precise  and  controlled  manner. 

The  bottom  most  plate  consists  of  a  irdwood  table  upon 
which  all  the  structure  above  is  either  directly  or  indirectly  fasten¬ 
ed.  This  bottom  plate  can  be  moved  vertically  by  means  of  four  Jack 
screws  whose  motion  ran  be  accomplished  by  means  cf  a  single  drive  shaft. 
This  drive  shaft  can  either  be  moved  manually  or  turned  by  a  motor.  The 
motor  employed  is  oi  the  gear  reduction  type  given  a  drive-shaft  speed 
of  5.7  rpm."  The  gearing  of  the  jack  screws  is  iuch  that  eight  turns 
of  the  motor  shaft  moves  the  table  one  inch.  The  vertical  position  of 
the  bottom  table  can  be  varied  over  a  range  of  some  eleven  inches. 

At  present  it  is  envisioned  that  the  vertical  degree  of  freedom 
will  be  used  to  position  the  various  probes  a*  the  correct  height  relative 
to  the  nozzle  axis.  The  radial  profiles  are  o  be  taken  in  the  horizontal 
direction.  If  it  proves  to  be  desireable  to  make  radial  profiles  using 
the  vertical  degree  of  freedom,  a  microswitch  has  been  positioned  relative 
to  the  drive-shaft  which  serves  the  function  of  automatically  stopping 
the  motor  after  a  half-turn  of  the  drive  shaft  or  after  any  multiple  there¬ 
of.  A  half-turn  of  the  drive  shaft  changes  the  vertical  position  of  the 
probe  0.0625  inches.  Dynamic  breaking  is  employed  so  that  there  is  no 
coasting  after  the  power  to  the  motor  is  interrupted. 

The  next  plate  up  is  made  of  alumin  *s  inch  thick  and  is 
constrained  to  move  in  the  axial  direction,  i.e.,  along  the  axial 
direction  of  the  Jet.  This  plate  rides  on  four  matching  precision 
roller  bearings  along  two  ground  case-hardened  steel  shafts  having 
a  nominal  diameter  of  1  inch.  The  tolerance  between  the  shaft  and  the 
bearing  is  0.001  inches.  This  plate  is  driven  by  a  25  rpm  gear-re¬ 
duction  motor  through  a  lead  Bcrew  having  a  -h  of  10  threads  per 
inch.  The  axial  motion  of  the  plate  can  be  d  over  a  range  of  32 

inches  so  that  without  moving  the  table,  an  profile  can  be  obtained 

which  starts  from  the  nozzle  exit  plane  and  ,ds  some  43  nozzle  dia¬ 
meters  (0.75  inch  diameter)  downstream.  The  *1  lead  screw  is  fitted 

with  a  microswitch  so  that  the  axial  motion  ct  be  made  in  units  of 
0.05  inches. 

The  next  plate,  also  made  of  aluminum,  is  3/8  inch  thich 
and  is  constrained  to  move  in  the  horizontal  di-ection.  The  slid¬ 
ing  apparatus  is  similar  to  that  for  the  table  below  and  differs 
c  „y  in  that  shafts  have  nomir.  1  diame  r  of  3/4  inch.  This 
plate  moves  relative  to  the  plate  below  by  means  of  a  hydraulic  piston 
which  has  an  8  inch  stroke.  The  purpose  of  this  piston  is  to  have 
the  capability  of  positioning  the  probe  and  the  jet  impignment  plates 
(described  below)  in  and  out  of  the  exhaust  jet  much  more  rapidly  than 
can  be  obtained  with  the  lead  screws  which  are  primarily  used  for  con¬ 
trolled  positioning.  A  displacement  transducer  is  mounted  in  such  a 
way  as  to  monitor  the  piston  travel  to  ensure  that  identically  the 
same  piston  stroke  occurs.  -32- 


Figure  21 — Probe  Positioning  Mechanism 


The  top  plate  upon  which  the  various  probes  are  to  be  fastened  Is 
made  of  aluminum,  3/8  inch  thich,  and  Is  also  constrained  to  move  in 
the  horizontal  direction.  The  sliding  arrangement  for  this  plate  is 
identical  as  for  the  plate  below.  This  plate  is  driven  with  a  14  rpro 
motor  and  a  lead  screw  having  a  pitch  of  20  threads  per  inch.  This 
lead  screw  is  also  fitted  with  a  microswitch  so  that  the  horizontal 
motion  can  b*  made  in  units  of  0.025  inches.  The  total  travel  is 
10.5  inches. 

A  control  panel  has  been  built  whichallows  the  probe  positioning 
table  to  be  operated  from  the  control  room  of  the  arc  facility.  Through 
the  use  of  counters  on  the  control  panel  the  position  of  the  probe  with¬ 
in  the  jet  can  be  established  from  the  control  room.  A  subsidiary  con¬ 
trol  panel  is  mounted  on  the  end  of  the  probe  positioning  table  to 
facilitate  probe  alignment  and  other  preparatory  procedures. 

B.  )  High  Pressure  Water  Supply 

A  high  pressure  gas-driven  water  supply  has  been  fabricated;  a 
schematic  is  shown  in  figure  22.  Basically  it  consists  of  a  pressure 
tank,  a  heat  exchanger,  a  collecting  tank  and  a  flow  meter.  The  system 
as  designed  will  provide  chemically  clean  water,  free  of  particles 
larger  than  5  microns,  at  operating  pressures  up  to  1000  psi.  The  pre¬ 
ssure  tank  holds  approximately  40  gallons  which  should  provide  a  run¬ 
ning  time  of  an  hour. 

C.  )  Mobile  instrumentation  Platform 

A  compact  mobile  instrumentation  platform  has  been  designed  and 
built.  All  auxilliary  equipment,  with  the  exception  of  the  probe 
positioning  and  coolant  systems,  required  for  operating  of  the  probes 
to  be  used. in  the  course  of  this  program  ic  provided  for  on  the  plat¬ 
form.  It  thus  provided  the  capability  of  moving  and  setting  up  a 
probe  experiment  in  the  shortest  possible  time.  The  instrumentation 
platform  is  shown  in  Figure  23.  Located  on  the  lower  shelf  can  be  seen 
the  constant  volume  tank,  the  vacuum  pump,  and  array  of  six  solenoid 
valves  (two  pair  of  which  are  arranged  in  series  opposition) ,  and  a 
thermos  bottle  which  provides  the  reference  ice  Junction  for  the 
thermcouple  which  reads  the  tomperature  of  the  gas  in  the  tank.  Located 
on  the  top  platform  are  gauges  to  read  probe  impact  pressure  and  vacuum 
tank  pressure  (to  be  replaced  by  transducers) ,  a  micrometer  needle  valve 
to  regulate  the  mass  flow  rate  of  gas  through  the  probe  and  into  the 
tank,  an  array  of  sample  bottles,  each  of  which  disconnects  for  analysis 
of  the  composition  of  the  gas  collected  therein,  and,  on  the  extreme 
right,  the  gas  inlet  to  which  the  probe  is  connected.  A  schematic 
diagram  of  the  lines  and  instrumentation  located  on  the  mobile  plat¬ 
form  is  shown  in  Figure  24.  Yet  to  be  added  to  the  assembly  are  a 
remote  switch  box  for  operation  of  the  solenoid  valves,  a  timer  con¬ 
nected  to  solenoid  #3,  and  an  inert  gas  bottle  for  flushing  probe  lines 
and  bringing  gas  sample  bottles  up  to  atmospheric  pressure  prior  to 
disconnect. 

In  operation,  vacuum  is  first  established  in  the  tank  and  the 
sample  bottle  (figure  24)  with  valves  one,  two,  five  and  six  open. 

The  initial  temperature  and  pressure  in  the  tank  are  recorded.  The 
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Figure  23 — Probe  Instrumentation  Platform 


valves  are  then  shut  and  the  probe  is  exposed  to  the  jet  to  be  in¬ 
vestigated.  The  impact  pressure  and  the  temperature  rise  of  the 
inner  coolant  with  no  gas  flow  (the  "tare")  are  recorded  at  this 
time.  Valve  three  is  then  opened,  allowing  a  sample  to  flow  through 
the  probe  (the  needle  valve  having  been  preset).  During  this  time, 
the  temperature  rise  of  the  inner  coolant,  the  temperature  of  the  gas 
at  the  probe  exit,  and  the  temperature  of  the  water  leaving  the  outer 
probe  cooling  jacket  (in  order  to  prevent  damage  to  the  probe) ,  are 
recorded.  At  the  conclusion  of  the  test,  valve  three  is  shut  and 
the  temperature  and  pressure  of  the  gas  collected  in  the  vacuum  tank 
are  noted.  Valves  one  and  two  are  opened,  thus  filling  the  sample 
bottle.  Valves  three  and  four  are  opened  to  purge  the  lines  with 
mertgas  and  to  bring  the  sample  bottle  to  atmospheric  pressure,  sub¬ 
sequent  to  which  the  bottle  is  disconnected  and  its  contents  analyzed 
for  the  ratio  of  0  to  N.  This  is  done  using  the  Abcor  Model  WB-7 
Radio  Frequency  Mass  Spectrometer  shown  in  Figure  25. 

The  enthalpy  of  the  gas  being  probed  is  found  from  the  following 
formulas 


Cpt  (  ATC  -  A\nc^) 


UO 


The  known  gas  composition  along  with  initial  and  final  values  of 
temperature  and  pressure  in  the  tank  and  the  time  of  the  run  allow  a 
simple  calculation  of  ftg  using  the  perfect  gas  law.  The  measured  gas 
exit  temperature,  gas  composition,  and  the  known  equation  of  state 
for  nitrogen-air  mixtures  provides  a  value  of  h2„.  The  other  quantities 
are  direct ly' measured  and  hlg  is  easily  determined. 


An  Avco  computer  program  for  the  thermochemical  properties  of 
nitrogen-air  mixtures  has  been  utilized  for  the  preparation  of  Figures 
26  through  30.  These  charts  given  dimensionless  enthalpy  vs  entropy, 
temperature,  density,  specific  heat  ratios,  and  speed  of  sound. 

Pressure  is  one  atmosphere  in  all  cases.  The  atom  ratio  of  oxygen  to 
nitrogen  is  a  parameter  on  all  curves.  Thus  the  known  values  of 
enthalpy  and  composition,  along  with  the  assumption  that  the  gas  m 
the  jet  is  in  equilibrium,  provide  values  of  temperature,  density , speed 
of  sound,  and  specific  heat  ratio.  The  measured  values  of  stagnation 
pressure  and  the  value  of  *  read  from  the  charts,  give  the  Mach  No.  from 
the  isentropic  flow  relationship. 


?* 

P 
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Finally,  the  Mach  No.  and  the  speed  of  sound  give  the 
The  assumption  that  the  static  pressure  in  the  jet  is 
to  one  atmosphere  will  be  checked  frequently  with  the 
probe  described  above. 


local  velocity, 
everywhere  equal 
static  pressure 
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Figure  25 — Radio  Frequency  Mass  Spectrometer 
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Figure  30--Nitrogen-Air  Mixtures!  Enthalpy  vs  Speed  of  Sound 
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D.  Pressure  and  Heat.  Flux  Distribution  Plates 


Plat  copper  plates,  nine  inches  in  diameter,  have  been  constructed. 
One  is  designed  to  determine  the  heat  flux  distribution  produced  by 
the  impinging  jet?  the  other,  the  pressure  distribution.  The  pressure 
distribution  plate  is  water-cooled  utilizing  a  number  of  concentric 
cooling  passages.  The  heat  flux  plate  is  uncooled,  and  is  instrumented 
with  a  number  of  transient  calorimeters.  The  pressure  distribution  plate 
has  undergone  preliminary  testing  in  the  Avco  Hydrogen -Oxygen  Rocket 
Motor  Facility  in  order  to  determine  its  ability  to  withstand  the  high 
heating  environment.  Pressure  taps  are  now  being  added  to  the  plate. 

A  total  of  eleven  taps  are  being  provided,  one  in  the  center,  seven  on 
one  radius,  and  three  others  to  insure  that  the  jet  and  plate  axes  of 
symmetry  are  coincedent . 

The  plate  designed  to  Jetermir.e  heat  flux  distribution  has  also 
been  completed.  Eleven  3/16‘'  transient,  copper  calorimeters  are  located 
in  such  a  manner  that  none  disturbs  the  radial  flow  passing  over  another. 
Each  calorimeter  consists  of  3/16  inch  diameter  and  l*s  inch  long  copper 
cylinder  into  which  a  thermocouple  has  been  brazed  m  place.  The 
thermocouple  consists  of  10  mil  chromel -alumel  wire  contained  in  a  2- 
hole  alumina  rod.  all  placed  in  a  stainless  steel  sheath.  The  thermo¬ 
couple  has  been  brazed  into  the  copper  slug  by  means  of  a  20  KW  vacuum 
induction  heater.  The  location  of  the  thermocouple  is  such  that  what¬ 
ever  temperature  the  thermocouple  senses  is  the  surface  temperature 
of  the  copper  slug  without  the  presence  of  the  thermocouple .  Once  the 
surface  temperature-time  history  is  determined  for  the  one  dimensional 
heating  of  the  copper  slug,  an  IBM  program  is  employed  to  determine 
the  heat  flux  to  the  copper  calorimeter  as  a  function  of  time.  This 
calorimeter  is  mounted  in  a  copper  sleeve  with  an  annular  air  gap  between 
the  calorimeter  and  the  sleeve. 

The  flat  plate  pressure  distribution  will  be  utilized  to  find  the 
radial  velocity  gradient  at  the  stagnation  point  of  the  plates,  Follow¬ 
ing  Bade<12),  the  velocity  distribution  produced  by  normal  impingement 
of  a  free  jet  upon  a  flat  plate  can  be  represented  approximately  by 
the  analytical  form 


u.  = 


w 


in  which  u.  is  the  component  of  the  velocity  parallel  to  the  surface 
at  the  edge  of  the  boundary  layer  r  is  the  radial  distance  from  the 
stagnation  point  (9  is  the  stagnation  point  velocity  gradient  and 
C‘i  is  a  coefficient  with  the  dimensions  of  an  inverse  length.  Using 
the  measured  pressure  data  and  assuming  isentropic  flow,  the  ratio 
of  local  velocity  to  free  stream  jet  velocity  {  U/U*.  ),  is  calculated 

for  each  port  location.  The  non-dimensional  parameter  K  (r)  *  ) 

\  u.  /  u. «.  )  is  computed  for  each  experimental  pressure  reading.  K(r) 

is  then  plotted  against  r  and  the  experimental  points  are  fitted  with  a 
straight  line  equation  of  this  line  is 

KCr)  =  (  /  u.«o)  (0  (  I  +  C,  r)  (5) 

•  h*)- 


The  intercept.  K(o)  of  the  fitted  line  at  r«o  therefore  gives  the 
coefficient  in  the  velocity  gradient,  relation 


(3  =  K(O)  {  U*o/T>j»0 


C« 


The  stagnation  point  velocity  gradient  so  determined  will  be  used  in  the 
Fay-Riddell (8)  equation  for  stagnation  point  heating  in  order  to  calculate 
the  heat  flux  to  the  plate.  Such  calculations  will  be  compared  with  the 
measured  values  of  heat  flux  obtained  as  described  above. 

III.  ANALYTICAL 

During  the  technical  information  meeting  of  21  July  1964, Aeronautical 
Research  Associates  of  Princeton  presented  an  analysis  of  turbulent  mix¬ 
ing  along  a  plane  free  jet  boundary.  Doth  incompressible  and  compressible 
cases  were  considered.  The  flow  model  and  nomenclature  used  in  the  analysis 
are  shown  in  Figure  31.  In  this  manner  it  was  possible  to  relate  all  of 
the  flow  properties  of  interest,  including  the  turbulent  exchange  coefficient, 
to  the  spreading  parameter  a*&/%,  a  quantity  which  can  be  experimentally 
determined.  An  important  assumption  made  in  this  analysis  is  that  the 
particular  form  of  the  assumed  velocity  profile  across  the  mixing  region 
does  not  greatly  affect  the  turbulent  exchange  coefficient  determination. 

Thus,  for  simplicity,  a  linear  profile,  U=  u6  +  3/S  was  chosen.  In  order 
to  provide  a  single  check  on  that  assumption,  an  alternate  sinusoidal 
velocity  profile  (which  should  be  chosen  to  reality),  was  chosen 
The  incompressible  results  obtained  under  each  assumption  are  shovJH  in  the  * 
table  beiow. 


VELOCITY  PROFILE 

Straight  Line 

Sinusoidal 

do/U* 

2/3 

.692 

*1. 

a/3 

3a/8 

n* 

-2a 

3 

-5a/8 

AJ./Us 

a/6 

a/8 

i* 

-.090a 

-a/16 

K 

.063a 

.067a 

It  can  be  seen  that  the  expression  for  K  is  changed  only  slightly 
by  the  alternate  velocity  profile  assumption,  thus  suggesting  that  the 
original  assumptions  of  a  straight  line  profile  is  sufficient  and  that 
K  is  not  sensitive  to  the  particular  profile  chosen. 

IV.  PLANS  FOR  THE  NEXT  QUARTERLY  PERIOD 

During  the  next  quarterly  reporting  period  it  is  intended  that  a 
good  portion  of  the  axial  and  radial  surveying  of  the  jet  be  completed. 


FtGoxE  3/  FLOW  MODEL  TOFJBt/L Ea/F  M/X/WG  /9a//9LYS/S 


This  will  be  preceded  by  completion  and  calibration  of  the  probea  to  be 
used,  and  final  construction  of  the  instrumentation  platform,  the  position¬ 
ing  mechanism,  and  the  high  pressure  water  supply.  Data  reduction  con¬ 
nected  with  these  surveys  should  also  be  substantially  completed.  It  is 
intended  that  a  computer  program  be  written  in  order  to  make  ouch  data 
reduction  fast  and  accurate,  requiring  only  directly  measured  quantities 
plus  values  read  from  the  nitrogen-air  mixture  diagrams  as  inputs. 
Aeronautical  Research  Associates  of  Princeton  will  be  consulted  regard¬ 
ing  the  calculation  of  turbulent  exchange  parameters  arising  from  these 
survey  results,  in  order  that  the  exchange  parameters  determined  may 
have  maximum  application  to  all  portions  of  the  program  and  to  future 
investigations. 

In  addition,  significant  progress  should  be  made  in  the  determin¬ 
ation  of  flat  plate  heat  flux  and  pressure  distributions.  These  quantities 
will  be  investigated  at  various  axial  positions  in  both  properly  ex¬ 
panded  and  underexpanded  Jets.  These  tests  will  be  preceded  by  completion 
of  the  revised  plates  and  further  testing  of  their  ability  to  withstand 
the  high  heating  environment  to  which  they  are  to  be  subjected. 
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SUMMARY 


Tlic  work  performed  by  ARAP  In  assisting  NRL  in  the 
monitoring  of  the  Acrothermal  Phase  is  reviewed.  Several 
experimental  studies  designed  to  provide  back-up  Infor¬ 
mation  for  the  more  extensive  efforts  of  the  other  groups 
Involved  are  described.  Recent  developments  In  ARAP's  own 
basic  research  program,  Including  the  completion  of  reports 
covering  the  stu.b-  of  jet  impingement  phenomena,  are  also 
discussed . 


UNCLASSIFIED 


LI 


..’.l.i 


.r  N'-jL‘Hor->  :c  *i-x  rji-3 


During  the  nix -month.  poxrl  od  just 
A  RAP's  over-all  program  oJ‘  research  .bed', 
of  aorot Dermal  processes  under  the  hhi! 
essentially  completed,  while  invests 
have  continued.  The  re  cult  a  of  an  ■.  ,.p.  ' 

impinging  jets  are  being  reported  1. 
iveport  i'foc,  j j,3  a; (j'l- ,  .no  ol  c ne n < 
ho will  be  distributed  in  the  near  > 


t: or  vice.  on 
,n  ■  o  aspects 
'.nan  in:.:, 
other  ?n  ear; 

stud;.'  of 
••;3  A Fi A:’ 
.iplete  ,  r.t/i 
ither  ore as 


.hi  v.'ii.u'h.  exporim.ui.cil  invoa  vignt.i  on  '  ■:<>.  .  :.:.uing  are 

t : .os t:  dealing  wit!,  v.he  phenomena  of  no'.....,  'lows  arid  fi 
problem  of  vehicle  hardening  bp  -none.,  o  wclor  for:;, 
pa.* ng,  'fhpnr.ef  1  cal  r hi. ini  c-f  of  f  •  '••«r  •  *  v‘v-  b 


extendeu  to  consideration  of  mix  I,.:;  of  yasoc-  in  or. 
effort  to  gain  a  better  under  eta..".  •>  '  ear  layer 

mixing  in  general  and  coupled  Hots  1  .  ;u  vdevilar, 
has  continued,  during  this  periou,  •  :  eio".  ;:f.L  in  th 
over-all  management  and  coordlnatio  .  o:  ;  . .  <cu":-:'.?rr.a 

f.'.ase.  As  a  part  of  thin  effort,  .-ml  ;;  .  .  tonal 

1 > '. ve s 1 1  >’;a 1 1  on s  Dave  been  carried  «. -v. n  p  .obi-  ■  n  of 
current  interest.  In  addition,  meet;!.:. ;  •  ...  :.:n.;  v. ;; 
participants  in  the  Aerothernal  fhase  La*-'  hoi..-.  :* 

order  to  set  program  goals  in  t:\. .  11,;'.  f.  ■>  •  •  .cat 

recant  developments . 


MANAGEMENT  AID’'  COOK.' 
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outgassing  of  the  bonding  material,  a  pressure  tap  and  a 
thermocouple  were  inserted  in  the  bonding  material  directly 
below  the  impingement  point.  Other  thermocouples  were 
provided  at  intervals  away  from  this  point  so  that  the 
temperature -time  history  could  be  recorded  throughout  the 
heated  region.  Earlier  tests  of  this  kind  which  were  not 
Instrumented  had  shown  that  a  failure  due  to  bond  decompo¬ 
sition  could  occur,  but  the  composite  material  in  that  case 
was  altered  to  Induce  such  an  effect  by  using  a  very  thin 
aluminum  layer.  On  applying  the  flame  to  the  instrumented 
sample  of  unaltered  material,  11  was  found  that  the  thermal 
deformations  which  were  apparently  of  different  magnitude 
in  the  aluminum  and  phenolic  nylon  layers  were  sufficient 
when  combined  with  the  decreased  bond  strength  nt  higher 
temperatures  to  cause  the  mechanical  separation  of  the 
layers  before  any  appreciable  pressure  build-up  could  occur. 

As  a  result,  a  further  test  was  planned  using  a  sample 
mounted  rigidly  in  a  frame  which  both  sealed  and  clamped 
the  edges  In  order  to  avoid  over-all  mechanical  deformation. 

It  is  felt  that  this  arrangement  will  simulate  better  the 
case  of  a  large  expanse  of  material  surrounding  a  spot  of 
localised  heating. 

Another  study  which  has  been  undertaken  under  the 
management  and  coordination  portion  of  ARAP's  effort  is  one 
designed  to  estimate  the  magnitude  of  certain  interference 
effects  which  may  be  encountered  in  the  course  of  the  free 
and  .impinging  arc -jet  studies  being  conducted  by  AVCO. 

Two  principal  effects  appear  to  be  of  possible  Importance. 

The  first  of  these  is  a  localised  effect  due  to  the  presence 
of  the  probe  itself  in  the  flow.  In  a  shear  flow  such  as  that 
In  a  jet,  a  probe  produces  local  flow  distortions  which  can 
result  in  erroneous  readings.  The  magnitude  of  such  errors, 
in  turn,  depends  upon  the  size  of  the  probe  relative  to  the 
flow.  Another  potential  source  of  interference  is  the  entire 
geometry  of  cue  test  setup.  In  tills  case,  the  whole  jet 
flow  field  can  be  distorted  as  it  passes  nearby  objects. 

This  effect  can  be  complicated  by  asymmetries  produced  in 
the  entrainment  flow  field  by  other  parts  of  the  apparatus. 

In  order  to  check  the  magnitude  of  such  effects,  ARAP  has 
begun  a  series  of  measurements  using  probes  of  the  AVCO 
design  in  several  sizes.  By  making  comparisons  between  each 
of  these  probes  and  a  reference  probe,  It  should  be  possible 
to  Isolate  the  probe  size  effect.  In  addition,  a  test  setup 
designed  to  simulate  AVCO's  will  be  used  and  compared  to  one 
of  different  geometry  In  an  attempt  to  evaluate  the  over-all 
Interference  effect. 
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During  the-  report  ip.,";  period,  three  meetings  took 
place  among  those  groups  engaged  in  the  Aerothermal  Phase 
(NHL.  AVCO,  OE,  and  ARAP).  Those  meetings  primarily 
involved  the  review  and  discussion  of  recent  results, 
current  activities,  and  immediate  future  goals.  The 
problem  of  evaluating  foam  packing  materials  as  hardeners, 
which  is  currently  under  study  at  UK,  was  discussed  at 
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Theoretical . 
of  free  shea 
period,  A 

moving  at  some  velocity  psstoc  over  a  cavity  of  stagnant 

jps  betv/if-n  fee  raovin, 
gas  causes  a  !:lx!  r.g  region  to  neve  lop  downstream  of  the 
initial  encounter  of  the  two  gases.  Expressions  have  been 
derived  and  programed  for  a  digital  computer,  and  results 
have  been  obtehogi  on  the  following  properties  of  such  a 
i  roe  shear  f±ow • 

1.  Tiie  1.0'' at the  dividing  streamline 

2.  Th*.  grow'.:*  of  f:.e  mixing  region 

g.  The  c:n  rg,  hro  inferred  a  ores;:  t'v  ..lining  region 

d ,  The  stagnation  pom.':  heat  *.!.,.;  cause  a  by  the  free 
oh  .  ur  flow  :i  tav.'i j.i.g  on  a  hi  mt  edged  wall . 


Initially,  this  t-erx  v.ar.  conf  j nod  to  mixing  of  a 
single  unif •. dubs.,  gv  fin.  mixing  of 

distinct  gasco,  or.*  in  tno*  main  stream  and  the  other  in 
the-  cavity  were  .incorporated  into  the  analysis.  The 
consideration  of  a  different  gas  .in  the  cavity  stems  from 
the*  possibi  11  of  decomposition  products  being  present 
duo  to  internal  heating.  Thus,  it  is  now  possible  to  solve 
for  the  abovo- ••i  icnt  lonf.d  characteristico  as  [‘vine  t  tons  of 
the  propcrt.'ir  s  o"  too*  two  gases  as  ve*ll  as  their  stagnation 
enthalpy  ratio  and  Much  number  of  the  main  flow. 

A  comprehensive  report  which  will  present  the  analysis 
and  results  for  a  number  of  oar.es  Is  in  preparation  and 
will  be  issued  as  a.  separate  document. 
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Experimental .  The  principal  part  of  ARAP's  experimental 
program  under  the  Aerothermal  Phase  lias  been  the  thorough 
investigal:  Lon  of  Jet  impingement  phenomena.  This  study  has 
been  particularly  Important  in  view  of  the  role  of  heat 
transfer  by  jet  impingement  for  both  small  hole  ( low  A/V  /'3) 
and  large  hole  or  "coupled"  (large  A/V^/3 )  conditions  of 
tiie  punctured  reentry  vehicle.  .During  the  past  six -month 
period,  the  primary  effort  has  been  directed  to  the  comp¬ 
letion  of  this  experimental  phase  and  the  preparation  of 
two  comprehensive  reports  in  which  the  detailed  results  are 
presented.  The  first  of  these  reports  is  complete  and  is 
to  be  Issued  as  ARAP  Report  No.  03.  It  covers  the  free  jet 
and  normal  impingement  studies  which  constitute  the  bulk:  of 
the  program.  A  second  report  covering  results  for  oblique 
impingement  is  to  be  issued  as  ARAP  Report  No.  64  and  is 
now  nearing  completion.  Most  of  the  pertinent  results  of 
these  studies  have  been  reviewed  in  previous  status  reports. 

In  view  of  the  increased  understanding  gained  in  recent 
months  of  the  general  nature  of  the  free  shear  layer  as 
applied  to  the  coupled  flow  problem,  it  has  been  decided  to 
place  less  emphasis  on  further  experimental  studies  in  this 
area  at  the  present  time.  However,  a  correspondingly  more 
thorough  investigation  of  the  behavior  of  foam  packing 
materials  has  been  planned  and  is  about  to  begin. 
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